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Preface. 



More than half a century ago Edouarb Roche wrote his celebrated paper on the 
form assumed by a liquid satellite when revolving, without relative motion, about a 
solid planet."^ In consequence of the singular modesty of Roche's style, and also 
because the publication was made at Montpellier, this paper seems to have remained 
almost unnoticed for many years, but it has ultimately attained its due position as a 
classical memoir. 

The laborious computations necessary for obtaining numerical results were carried 
out, partly at least, by graphical methods. Verification of the calculations, which as 
far as I know have never been repeated, forms part of the work of the present paper. 
The distance from a spherical planet which has been called ^^ Roche's limit'' is 
expressed by the number of planetary radii in the radius vector of the nearest 
possible infinitesimal liquid satellite, of the same density as the planet^ revolving so 
as always to present the same aspect to the planet. Our moon, if it were homo- 
geneous, would have the form of one of Roche's ellipsoids ; but its present radius 
vector is of course far greater than the limit. Roche assigned to the limit in 
question the numerical value 2*44 ; in the present paper I show that the true value 
is 2*455, and the closeness of the agreement with the previously accepted value 
affords a remarkable testimony to the accuracy with which he must have drawn his 
figures. 

He made no attempt to obtain numerical solutions except in the case of the 
infinitely small satellite. In this case the figure is rigorously ellipsoidal, but for 
a finite satellite this is no longer the case ; nor do his equations afford the means of 
determining exactly the ellipsoid which most nearly represents the truth. These 
deficiencies are made good below, and we find that even in the extreme case of two 
equal masses in limiting stability the ellipsoid is a much closer approximation to 
accuracy than might have been expected. 

It is natural that Roche, writing as he did half a century ago, should not have 
been in a position to discuss the stability of his solutions with completeness, and 
although he did much in that direction he necessarily left a good deal unsettled. 

In 1887 I attempted the discussion of some of the problems to which this paper is 
devoted, by means of spherical harmonic analysis, t Poinoare's celebrated memoir 
on figures of equilibrium^ was published just when my work was finished, and I kept 
my paper back for a year in order to apply to my solutions the principles of stability 

■^^ ^*La figure d'une masse fluide soumise a rattraction d'un point eloigne," 'Acad, des Sci. de 
Montpellier,' vol. 1, 1847-50, p. 243. 

t ''Figures of Equilibrium of Rotating Masses of Fluid," 'Phil. Trans. Eoy. Soc.,' vol. 178 (1887) A, 
pp. 379-428. 

J "Sur Pequilibre d'une masse fluide animee d'un mouvement de rotation," ' Acta. Matli.' 7 : 3, 4 (1885), 
pp. 259-380. 
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enounced by him. The attempt is given in an appendix to my paper, but 
unfortunately I failed to understand his work completely, and my investigation, 
as it stands, is erroneous from the fact that one term in the energy is omitted. "^ 
I think, however, that the defect may easily be made good. 

The analysis of the present paper is carried out by means of ellipsoidal harmonic 
analysis. In the course of the work it becomes necessary to refer to previous papers 
by myselfj, all published in the 'Philosophical Transactions'; they are : '^ Ellipsoidal 
Harmonic Analysis," vol. 197 (1901) A, pp. 461-557; " The Pear-Shaped Figure of 
Equilibrium of a Rotating Mass of Liquid," vol. 198 (1902) A, pp. 301-331 ; ''The 
Stability of the Pear-Shaped Figure of Equilibrium, &c.," vol. 200 (1903) A, 
pp. 251-314; " The . Integrals of the Squares of Ellipsoidal Surface Harmonic 
Functions,^^ vol. 203 (1904) A, pp. 111-137. These papers are hereafter referred to 
by the abridged titles " Harmonics,'' " The Pear-Shaped Figure," " Stability,'' and 
" Integrals." 

The analysis involved in the investigation is unfortunately long and complicated, 
but the subject itself is not an easy one, and the complication was perhaps 
unavoidable. 

The principal inducement to attack this problem was the hope that it might throw 
further light on the form of the pear-shaped figure in an advanced stage of develop- 
ment when it might be supposed to consist of two bulbs of liquid joined by a very 
thin neck. The arguments adduced below seem to show that such a figure must be 
unstable. 

M. LiAPOUNOFF has recently published a paper in which he states that he is able 
to prove the instability of the pear-shaped figure even when only infinitesimally 
furrowed.*!* In view of my previous work on the stability of this figure, and from 
other considerations it seems very difficult to accept the correctness of this result. 

At the end a summary is given of the conclusions arrived at, and this last subject 
is discussed amongst others. 



Part I. — Analysis. 

§ 1. The Stability of Liquid Satellites. 

This paper deals with two problems concerning liquid satellites which possess so 
much resemblance that I did not for some time perceive that there is an essential 
difference between them. One of these is the determination of the figures and of the 
secular stability of two masses of liquid revolving about one another in a circular 
orbit without relative motion of their parts. We may refer to this as the problem of 

^ This is the term denoted Jw^ (8iy/I below. 

t "Sur un probleme de Tchebychef," 'Acad. Imp. des Sci. de St. Petersbourg,' vol. 17, No. 3 (1905). 

Y 2 
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'Hhe figures of equilibrium"; the other may be called '^Roche's problem/' and it 
differs only from the former in that one of the two masses of liquid is replaced by a 
particle or by a rigid sphere. However, in the numerical solutions found hereafter, 
Roche's problem is slightly modified, for the rigid sphere is replaced by a rigid 
ellipsoid of exactly the same form as that assumed by the other mass of liquid in the 
problem of the figures of equilibrium. Thus, with this modification, the two problems 
become identical as regards the shape of the figures ; but, as we shall see, they differ 
widely as to the conditions of secular stability. This difference arises from the fact 
that in the one case there are two bodies which may be subject to tidal friction, and 
in the other there is only one. 

If in either problem there is no solution when the angular momentum has less than 
a certain critical value, if for that value there is one solution and for greater values 
there are two, then the principle of Poincare shows that the single solution is the 
starting point of a pair of which one has one fewer degrees of instability than the 
other. If, then, one of the two solutions is continuous with a solution which is clearly 
stable, it follows that the determination of minimum angular momentum will give us 
the limiting stability of that solution ; and this is the point of greatest interest in all 
such problems. 

Our two problems differ in the value of the angular momentum of which the 
minimum has to be found. For, if in Roche's problem the second body is a particle, 
it has only orbital momentum ; if the second body is a sphere, it must be deemed to 
have no rotation ; and, finally, in the modified form of the problem, the rotational 
momentum of the rigid body must be omitted from the angular momentum, which 
has to be a minimum for limiting stability. 

It will be useful to make a rough preliminary investigation of the regions in which 
we shall have to look for cases of limiting stability in the two problems. For this 
purpose I consider the case of two spheres as the analogue of the problem of the 
figure of equilibrium, and the case of a sphere and a particle as the analogue of 
Roche's problem. 

Let p be density, and let the mass of the whole system be |7rpa^ ; let the masses ot 
the two spheres be f7rpa^X/(l+X) and ^Trpd?l(l^\), or for Roche's problem let the 
latter be the mass of the particle. 

Let r be the distance from the centre of one sphere to that of the other, or to the 
particle, as the case may be ; and co the orbital angular velocity. 

In both cases we have 

The centre of inertia of the two masses is distant r/{l+X) and Xr/(lH-X) from their 
respective centres, and we easily find the orbital momentum to be 

X 



^7rpa>^(or^ 



(1 + X)^' 
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In both problems the rotational momentum of the first sphere is 

|(|.,a>^(^^J - 

In the first problem the rotational momentum of the second sphere is 

and in the second problem it is nil. 

If, then, we write L^ for the total angular momentum of the two spheres, and Z 
for that of the sphere and particle, we have 



2 



JL/i 






4 ^5 



^ITpSi^O) 



l+\'^' 






\r' 



U{l + Xf' (l + X)^a^J 



^ 



6/3 






\r' 



L^(l + X)^/« (l + X)^a^J 



J 



• • 



* * 



(1) 



On substituting for a its value in^ terms of r, these expressions become 



/siY' 



A = ii^pr-^. [I (1 +^^'^) (1 +^)^'^ (^) 



r 



^/^ r +X - 



\v^; 



^/ J 



JL/c 



{i^pfri 



a^ 



(1+xy 



a 



3/2 



r 



\l/2" 



ix-(i+xr(| +x L^ 



To determine the minima of these functions, we differentiate with respect to r, and 
equate to zero. 

Then, if ri, rg denote the two solutions, we find 



7\ 

a 
a 



6_ 
5X 



(1+X^/^)(1 + X)^/^ 



= fx^/^(i+xp. 



Whence 



/4 



Minimum A = (|^pf^a^ 4 {-^jf' (i^tp^ ' 



Minimum 



X7/6 



r — /^iTrnW^a^ 4 /- 2 \i/4 ___^___ 
^2 -- V3^P; «*-'*=U35/ d + xV^^^^' 



. . / 1 \l/4 

The ratio of the minimum of A to that of L2 is (-^+1 ] . Thus as X rises from 

to 00 this ratio falls from infinity to unity. 

All the possible cases of the first problem are comprised between X = and X = 1. 
When X = 0, ri = oo; and when X = 1, 



r 



a 



1 = y(it,2^/^) = 1-738. 
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Thus, in the problem of the figures of equiUbrium, if one of the two masses is large 
compared with the other, the two must be far apart to secure secular stability. This 
is exactly what is to be expected from the theory of tidal friction, for limiting 
stability is reached when there is coalescence of the two solutions which correspond 
to the cases where each body always presents the same face to the other.^ 

The result when the two masses are equal becomes more easily intelligible when it 

is expressed in terms of the radius of either of them. That radius is a/2^'^^, so that 

when X = 1 

' a 



r = l-738a = 2*191 



21/3 



Thus, in the latter case, limiting stability is reached when the two spheres are 
nearly in contact with one another, for if r were equal to twice the radius of either 
they would be touching. 

When the two bodies are far apart, the solution may be obtained by spherical 
harmonic analysis, and has comparatively little interest. But when the bodies are 
equal or nearly equal in mass, limiting stability for the figure of equilibrium would 
seem, from this preliminary investigation, to occur when they are quite close together. 
Accordingly, in finding numerical solutions hereafter, I have devoted more attention 
to this case than to any other. 

Turning now to the solution of the analogue of IIoche's problem, we see that when 
X = 0, ^2 ^ 0. This would mean that a very small liquid satellite could be brought 
quite up to its planet without becoming unstable. But we shall see that, when the 
satellite is no longer constrainedly a sphere, instability first occurs through the 
variations in the shape of the satellite. This preliminary solution does not, therefore, 
throw much light on the matter, excepting as indicating that we must consider the 
cases where the satellite is as near to the planet as possible. 

Next, when X = 1, we have 



^2 



y(f.2V3)a=yi,^.(^)=. 1-549 



a 



Thus, when the two masses are equal^ their distance apart is only about 1^ radii of 
either, and they will overlap. Here again it would seem as if stability would persist 
up to contact, but, as before, instability first sets in through variations in the shape 
of the satellite. 

Finally, when X is large, rg also becomes large. This case is the same in principle 
as that considered in the pi^oblem of the figures of equilibrium, for it means that if a 
large liquid body (formerly called the satellite) be attended by a small rigid body 
(formerly called the planet), secular stability will be attained when the small rigid 
body has been repelled by tidal friction to a great distance from the large liquid body. 

* See 'Eoy. Soc. Proc,,' vol. 29, 1879, p. 168, or Appendix G (b) to Thomson and Tait's 'Nat. Phil' 
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As this case may be adequately treated by spherical harmonic analysis, it need not 
detain us, and we see that the most interesting cases of Roche's problem are those 
where X lies between and 1. 

§ 2. Figures of Equilibrium of a Rotating Mass of Liquid and their 

Stability. 

A mass of liquid, consisting of either one or more portions, is rotating, without 
relative motion of its parts, about an axis through its centre of inertia with angular 
velocity o). We choose as an arbitrary standard figure one which does not difier 
very widely from a figure of equilibrium, and we suppose that any departure from the 
standard figure may be defined by two parameters e and /, which may be called 
ellipticities. It is unnecessary to introduce more than two ellipticities, because the 
result for any number becomes obvious from the case of two. We also assume a 
definite angular velocity for the standard configuration. 

Let F{e, /) denote the gravitational energy lost in the concentration of the system 
from a condition of infinite dispersion into the configuration denoted by e, f 

Let I{e^f),a){e^f) denote the moment of inertia and angular velocity about the 
axis of rotation in the same configuration. 

The initial values of these quantities are those for which e =/= 0, and are 
F(0, 0), J (0,0), CO (0,0). These all refer to the arbitrary standard configuration; 
they are therefore constants, and I shall write them F, /, m for brevity. 

Let ellipticities e, / be imparted to the system, and let the angular velocity be so 

changed that the angular momentum remains constant. 

Then 

I(e,f)a){e,f) = 1(0, 0) co (0, 0) = Ico, 

The kinetic energy of the system is half the square of the angular momentum 
divided by the moment of inertia ; and since the angular momentum is constant it 
is equal to ^ {Io)fjI (e, /). 

Thus the whole energy of the system, both kinetic and potential, is equal to 



If F(e,/) = F+SF, I{e,f) = I~{-SI^ the expression for the energy as far as 
squares of small quantities is 

~-(F+8F) + iil|i' = --V+Uco'--SV-~'W8I+W^Z 

The first two terms may be omitted as being constant and of no interest, and the 
energy with the sign changed, so that it is the lost energy of the system, becomes 
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Since a; is constant, we may write this 



S(F+|a>^/)~Ja)^(^'. 



On developing this by Taylor's theorem, it becomes 



4>/4+i^^li+^/^^+i/^l^j(^+i-^^^^ 



o)V 3/ . .dP'^ 






The condition for a figure of equiUbrium is that the first differentials of the energy 
with respect to the ellipticities shall vanish. If, therefore, eo, /o denote the equilibrium 
ellipticities, the equations for finding them are 



'^2 



\2 \ 



(O 



de 3e 



,^,^;{VH^'i)- , 






de 



aV^W^-^"f''<''^*"'^> 



CO 

7 






= 0, 



0. 



Multiplying the first of these by e and the second by/, and adding them together, 
we find 



el+f^){V+WI) 



de 



^2 C^2 Ci2 

^^0J^2 + (^/0 + 60/)^g^.+//0g^ 



(v+wi) 



+ 



OJ 



2 



I 



eeo ( — ) + (e/o + eo/) YeTf ■^j\df) 



de 



On substituting this in the expression for the lost energy, it becomes 



92 ga g2 

1 e (e - 2eo) -^ + (e/- e/ - Co/) to. + i/ (/- 2/) ^^^ 



(F+|a.^i) 



2 



/ 



. {e-2eo) ('1^ J+ 2 {ef~.ef,-ej)f^ |. +/(/-2/o) (|.J" 



Now let §6, S/ be the excesses of e and jT above their equilibrium values 60, ^, so 
that e = eo + S^, y = yo + ^ Then on substitution in the expression for the lost energy 
it becomes 



1 {(8e)^_e„^}i_ + {Sea/-6„/o} -,-|-^+ i mr-f;'} 






dedf 



df] 



{V+^c^'I) 



0) 

7 



'{{Sey-e,'} (l^y + 2 {Se8/-eo/o} |^ f^ + {(8/^-/0^} (|.)" 



Since 60, ^o are constants, the portion of this involving 60, fo explicitly is constant, 
and may be dropped. 
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Thus the variable part of the lost energy may be written 



1 (Sey 



2 






de' 



2 



de' I\de 



+ 868/ 



d'V 
Idedf 



+ W 



av 2 di di\ 



dedf I de df/ 



+ J (¥)' 



Lgy. + 2«^3y. j^l^gy^ 



\2^ 



This is a quadratic function of the departures of the ellipticities from their 
equiHbrium values, and the form is obvious which the result would have if there 
were any number of ellipticities. 

Since the condition for secular stability is that the energy shall be a minimum, 
the lost energy must be a maximum, and therefore this quadratic function of 
Se, 8/, &c. , must always be negative in order that the system may possess secular 
stability/^ 

If i^ is a quadratic function of n variables, Xi^ x^, x^, &c., so that 

-r ^^22*^2 "•" '"^23'^2'^3 ' • • • 

+ 0^33.^3 + . . . , 

it is known that the condition that F shall always be negative for all values of the 
variables is that the series of functions 



a 



11 ? 



a 



11? 



a 



12 



^12? ^22 



a 



11? 



a 



12? 



a 



12? 



a 



22y 



a 



13 



a 



23 



^13? ^23? ^33 



• 3 



• (2) 



shall be alternatively negative and positive. 

Since we might equally well begin with any one of the variables, it follows that 

<^ii? ^22^ (^S8'" must all be negative; also ^12^—^11^22, ctii — <^ii<^^z> %3^—<^22<^33--- must all 
be negative if F is always to be negative. 

Now, suppose that F is the function of lost energy for a system with ^ + 1 degrees 
of freedom, but that a constraint destroys one of the degrees. If the system has 
secular stability, the n determinants must have their appropriate signs, and when the 
constraint is removed, the new additional determinant must have its proper sign in 
order to secure secular stability. It follows that stability can never be restored by 
the removal of a constraint if the system was unstable when the constraint existed; 
but stability may be destroyed by the removal of a constraint. 

* This result is also given, but with less detail, in my paper on "Maolaurin's Spheroid," in 'Trans. 
Amer. Math. Soc.,' vol. 4, No. 2, pp. 113-133 (1903). 
VOL. OOVL— A. Z 



170 



Sm G. H. DAEWIN ON THE 



§ 3. On the Possibility of joining Two Masses of Liquid by a thin Neck 

This whole investigation was undertaken principally in the hope that it might lead 
to an approximation to the form of the pear-shaped figure of equilibrium of a rotating 
mass of liquid at the stage when it should resemble an hour-glass with a thin neck. 
It seemed probable that such an approximation might be obtained in the following 
manner : — 

Two masses of liquid are revolving in an orbit about one another without relative 
motion of their parts, so that they form a figure of equilibrium. Imagine them to be 
joined by a pipe without weight, through which liquid may flow from one part 
to the other. A flow of liquid will in general take place between the two parts, but 
there should be some definite partition of masses, corresponding to a given distance 
apart, at which flow will cease.. At this stage we should have an approximation 
to the hour-glass figure of equilibrium. 

In this section a special case of this problem is considered, in which the detached 
masses, to be joined by a pipe, are constrained to be spheres. 

If the notation of § 1 be adopted, it is clear that the system is defined by the two 
parameters r and X. In accordance with the notation of § 2 we denote the lost energy 
of the system by V and the moment of inertia by /. It is easily shown that 



V = (|^/>)V 



_X__ a 3 Jj-X^ 
l{l+Xf r^{l+Xf'\ 



For brevity write 



I 



4_„„5 

■gTToa 



r 






i2 ^2 ' 5 



F 



k 



{i+xy 



G 



(l+Xf" 



and let F\ G', F" , G" denote their first and second differentials with respect to X. 
The equations for determining the configuration of equilibrium are 



37 , „a/ 



dr 



dr 



0. -^r—+^(0 



ax 



3X 



0. 



The first of these gives at once 



2 4 / ^ \ 



For determining the form of the second we have 



V = (|,rp)V 



F^+iG 

7' 



I " -fTTpa' 



a^ * 
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If we differentiate these with respect to X, substitute in the second equation of 
equilibrium and give to co^ its value in terms of r^, we find that the result is 



a' 



+ 



15 ±_ ^ 

' G'r 



+ 3 = 0. 



Now 



F' 



l-\ 

{1+xr 



G' 



5 K'l^-l 



and 



G 



3. 
5 






Hence the equation for determining r for a given value of X is 



a^«_9_JL_+XV3+x2^ a 
r" '(l + X)V«(l+Xi/«)r 



+ 3 = 0. 



(3). 



This cubic has 'three real roots of which one is negative and has no physical 
meaning ; the second gives so small a value to r that the smaller sphere is either 
wholly or partially inside the larger one. The third root is the one required. 

In order to present the result in an easily intelligible form it may be well to express 
it also in terms of the radius of the larger of the two spheres, say ai, where 



a 



3 



a' 



l+A 



The following is a table of solutions for various values of X : 



X^i\ 


r/a. 


r/ai. 


r/ai - (1 + X^i% 


0-0 


1-304 


1 • 304 


0-304 


0' 


1 




323 




323 





•223 


0' 


2 




368 




371 


0- 


171 


0' 


3 




426 




438 





•138 





4 




486 




"517 





"117 





"5 




543 




■604 


0- 


•104 





6 




590 




697 





"097 





•7 




"625 




"793 





093 





"8 




"649 




•893 





093 





"9 




•662 




•995 





"095 


1-0 


1-666 


2-099 


0-099 



The solution is exhibited in fig. 1, the larger sphere being kept of constant size and 
the successive smaller circles representing the smaller sphere. Many of the circles 
pass nearly through one point, and it has not been possible to complete them without 
producing confusion. 

The fourth column of the table gives the excess of r above the sum of the two radii 
of the spheres, and it shows what interval of space is unoccupied by matter. It is 
remarkable how nearly constant that interval is throughout a large range in the 
values of X. 
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Fig. 1. Solutions for two spheres of liquid joined by a weightless pipe, for successive values of X^^^. 

In the case where the two bodies are no longer spheres, the equation corresponding 
to the cubic (3) becomes very compHcated. It is therefore desirable to discover 
whether in any mven solution of the figure of equilibrium the two detached masses 
are too far apart to admit of their being joined by a weightless pipe, or whether they 
are too near. This may be discovered in the following way : — 

Let Tq be the solution of/(a/ro) = 0, where 



l+X^^'+X^/^ a^^ 

2-TZ TTTJ^-TT: TTTTT - + O. 



/» / a \ a 9 

•^ \r/ r"' ~ ^ (1 + Xy'' (1 + V'«) r 



(4). 



There is only one solution of /" = between r equal to infinity and the case when 
the two spheres touch. Hence we can determine on which side of Tq any given value 
of r lies by merely considering whether f changes from positive to negative or from 
negative to positive as r increases through the value Tq. 



a 



a 



a 



Now if -^ + S ( ^ I be any neighbouring value of — , we have approximately 

^0 \^/ '^ 



/ 



a 



= 3 



a 



JL "J" A. "1" A 



in' '{l + Xf''{l + X'^')^ \r 



81- 



If we express aY^o* in terms of ro/a by means of the equation for r^, this may be 
written 



•^[rj (l + Xf" 



1 + X'" 



nn 



XV3 



a 



ai 1 + X^''' 



S -, 



r 



where as before 



ai 



a' 



l + X* 



Now the fourth column of the table shows that l+X^/*^— r/a^ is negative, and the 
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last term inside the bracket is also negative. Hence if S ( - ) is negative, fi-) is 



a 



a 



positive, and vice versd. But if S ( - is negative, r is greater than Tq and the two 



T 



masses are too far apart to admit of junction, and vice versd. 



a 



Therefore if for a given solution for detached masses f\-) is positive, the masses 



r 



are too far apart to admit of junction by a weightless pipe, and if it is negative they 
are too near. 

When in the general case we form a function y(a/r), such that when the ellipticities 
of the two masses are annulled, it reduces to the above function, its sign wdll afford 
the criterion as to whether the masses are too far or too near to admit of junction by 
a thin neck of liquid. I return to this subject below in § 13. 

The solution of the problem when the two masses are constrainedly spheres is so 
curious that it seems worth while to consider its stability. This may be done by the 
method of § 2. 

The system depends on two parameters r and X, and the stability will depend on 
three functions, which are defined as follows : — 






2 /9J\2 



0- — 

9r* I \dr 



^'^'^ drdX """^ drd\ IdrdX' 



dw 



du 



(O 



2 /9JN2 



{X,\} = ^ + W^,- 



These functions correspond to ^n, c^i2j <^22 of (2) in § 2, and we see that for secular 
stability {r, r} and {X, X,} must be negative, and 

A = {r,r}{X,\}-l{r,X}J 
must be positive. 

Without giving the details of the several differentiations, I may state that if we 

write 



H = (Ittp) V . 



q3 



.2 



F-„ + m 



a 



so that H is essentially positive, we find 






H 



-2r(y^'^, + f(?'); 



H 



F 



r 
a' 



iF"-- + lG"(l-^d 



a' 



a^ ^ 



a^ ^ 



>• 
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From the equation 






we find, after some reductions, 



On substituting for F\ G' their values in terms of X, I find that this expression 
reduces to 

X^/^(1+X)^a^' 

an essentially positive quantity. 
On substitution in A I find 



A=^(^^:+i6' 



!)fe^'^(*''-^?)-K^'"5+*'^'J" 



(5) 



The factors outside [ ] are essentially positive and do not aflfect the sign of A, and 



n/.2 



it is clear that A can only be positive if |- (r— i^-^ is positive. But A must be positive 

a 



n/.2 



for secular stability; hence stability can only be secured by jG—F—^ being positive, 

a 

and it is not necessarily so secured. But if this function is positive, so also is {?', r}, 

and if this last is positive the system is unstable. Hence stability is always 

impossible. As a fact, in all the solutions given above {r, r} is positive, and we 

should have to move the spheres much further apart to make it negative, and 

therefore on this ground alone the system, is always unstable. But A is sometimes 

positive and sometimes negative and vanishes for a certain value of X. As the 

vanishing of A puzzled me a good deal, I propose to examine the matter further. 

Before doing so, however, I will show that the instability of the system may be 
concluded from other considerations. 

It was proved in § 1 that two spheres, unconnected by a pipe, are in limiting 
stability when their distance apart is given by 



r^ 6 /, . xK/.'^X/. . xV1/.q A (~ 






&' 



= ^{l + X«'«){l + X)V3 = |^|,. 



This is the condition that {r, r} should vanish. 

When X is zero the two spheres in limiting stability are infinitely far apart, and 
when X is unity they are as near as possible, and r = l'738a. 

Now the table of solutions in the case where the two are connected by a pipe shows 
that they are furthest apart when X is unity, and that then r = l'666a. 

The removal of the constraint of one degree of freedom may destroy stability, but 
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cannot create it. Hence, when two spheres revolve about one another, the opening 
of a channel of communication between them may destroy stability, but cannot create 
it. When two equal spheres revolve about one another at such a distance that they 
could be connected by a pipe and yet remain in equilibrium, their distance is 1*666 ; 
but they are then unstable, because 1*666 is less than 1*738. The opening of a pipe 
between them, being the removal of a constraint, cannot make the motion stable. 
A fortiori the like is true when the two spheres are unequal in mass. 

Hence the system of equilibrium of two spheres joined by a pipe is unstable in all 
cases, 

I will now consider the meaning of the vanishing of A. 

Having evaluated the angular momentum of the system corresponding to the 
several solutions tabulated above, I found it had a minimum when X^^^ = 0*254. 
Such a solution is a critical one and is the starting point of two solutions of which 
one must have one fewer degrees of instability than the other. The vanishing of A 
must have the same meaning, but it remains to be proved that minimum angular 
momentum is secured by the vanishing of A. 

The angular momentum is /w, and is therefore proportional to /x, where 

On equating -J^ to zero so as to find its minimum, we have 



r^ . o ^/\ . / TTT r^ a^\ dr 



a^ 



2 ' 5- / M- g^2 ^-JdX 



Now since 



we have 



^*' _^ i _!_ 5 ^'' ^' _ A 
^+3 + 2^,^-^. 



dr _ >-ir(l-V/y X~^/^ (1+X) 



-4 



dX 



^ a 



On substituting this value in the equation dix/dk = 0, I find that the result may 
be written 

X^/^(l + X)^an' aV^ aH a^^ ^^ A' a + ^ / " ^- 

The first term of this is the same as the first term inside the bracket in the 
expression for A in (5). On comparing the two second terms together we see that 
A = is the condition for minimum angular momentum, if 
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that Is to say, if 



or If 



^^(*''i+^')-^^' ?+««'=»• 






But this last Is true, being the equation (3) determining the figure of equIHbrlum ; 
hence A = gives minimum angular momentum. 

Since two liquid spheres cannot be joined stably by a pipe, it seems very improbable 
that two tidal ellipsoids could be so joined as to become stable. Indeed, if the 
distortion of the surfaces of the two masses Into ellipsoidal forms may be regarded as 
due to the removal of constraints whereby they were previously maintained In a 
spherical form, stability Is Impossible. 

The question as to whether or not there is an unstable figure with a thin neck will 
be considered later, for the present we are only concerned with the conclusion that 
there Is no stable figure of this kind. 

Mr. Jeans has treated an analogous problem In his paper on the equilibrium of 
rotating liquid cylinders,"^ and has concluded that the cylinder will divide stably Into 
two portions. The analogy is so close between his problem and the three-dimensional 
case, that It might have been expected that the analogy would subsist throughout ; 
nevertheless, If we are both correct there must be a divergence between them at some 
point. 

§ 4. Notation. 

As the solution given below Is effected by means of ellipsoidal harmonic analysis, 
it Is well to state the notation employed. It Is that used In four previous papers to 
which references are given In the Preface. 

In ''Harmonics" the squares of the semi-axes of the ellipsoid were 



a" = fIv'-^ Y^j , h' = ¥ (v'^ 1), c' = kV, 



The rectangular co-ordinates were connected with ellipsoidal co-ordinates v, /x, <^ by 

Z^ 2 2 1— iScos2(i 
* 'Pliil. Trans. Eoy. Soc.,' Series A, vol. 200, pp. 67-104. 
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The three roots of the cubic 

/y»-s ~.a ^a 



+ T#-+-.f— == 1 



a^-\u ¥-\-u c^-i-u 
were 

1— p 

Lastly p ranges from oo to 0, /x between ± 1, <^ from to 2tt. 
In the two later papers, I put 

K^ = T^ , K^ = 1—K^, V = — ; — , iut = sin ^ ; 

and for convenience I introduced an auxiliary constant /3 (easily distinguishable from 
the ^ of the previous notation) defined by sin ^ = /csiny. 
The squares of the semi-axes of the ellipsoid were then 

a^ = ^' ^Q^' y jy2 ^ /g' cos^ /3 ^2 ^ /^' ^ 
sin^ ft ' sin^ ft ' sin^ ft ' 

The rectangular co-ordinates became 

^ = P ^ ^ ( 1 — K^ sin^ ^) cos^ (f), '4 = -?~o-^ cos^ ^ sin^ S, ^x = -r-^-?, sin^ 6(1— k^ cos^ ^X 
F sm^yS ^ F sm^^ ^ F snr^ , 

The roots of the cubic were 



sin^ /3 K" 



tfi = . g 5 U2 = k^ sin^ ^5 % = --^ (1 — K ^ cos^ <jf>). 



The notation employed for the harmonic functions is that defined in '' Harmonics." 

§ 5. The Determination of Gravity on Roche's Ellipsoid. 

In Roche's problem a mass of liquid, which assumes approximately the form of an 
ellipsoid, revolves in a circular orbit about a distant centre of force without any 
relative motion. In the present section it is proposed to evaluate gravity on the 
surface of this ellipsoid. I intend to solve the problems of the present paper by 
means of the principles of energy, and for that purpose it is necessary to determine 
the law of gravity. 

Suppose that the ellipsoid of reference, defined by Vq, is deformed by a normal 
displacement defined by the function pf ((jl, <^), where p is the perpendicular from the 
centre on to the tangent plane at /x, <^. This deformation must be expressible by a 
series of ellipsoidal harmonic functions, and therefore we may assume 

VOL. CCVl.^ — A, 2 A 
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The typical term written down must be deemed to include sine-functions as well as 
cosine-functions, and all those types which I have denoted by P, C? S in ^' Harmonics." 

On multiplying each side of our equation by any harmonic function, and integrating 
over the surface of the ellipsoid, an element of which surface is denoted by dcr^ we 
find in the usual way 

/(/!-, jy) P/ (/x) C/ {<j))pd(r 



p. = 



m{f^)^n^)jpd<r 



Suppose that / (jn, <^) is zero everywhere except over a small area Ba situated at 
the point ix\ ^^, and that it is there equal to a constant c ; also let p^ be the value 
of p at this area Scr, 

Then the mass of the inequality is 

where p is the density of the solid ellipsoid which is deformed. 

Next let us suppose that the mass of the inequality is unity, so that 

cp^pSa = 1. 
Then we have 

f/(i^. </>) f / (i^) (^i {<i>)pdcT - cf / (/) c/ (f )/Sa = If/ (/.o c/ (f )• 

J p 



Hence 



e/ = 



1/ (/^O ^i W) 



p 



[f/WC/(<^)Ti9(^0- 



I now write M for the mass of the eUipsoid, and shall subsequently make it equal 

to f Trpa^ - — r- , while the mass of the distant particle will be -^ or ^irpd? — — . 
^ 1+X ^ X ^ ^ 1-i-X 



Since 



ppda- = 3lf, and f o {p) Co (^) = 1, we have e^ = ~-^. 



Thus an inequality representing a particle of unit mass at /x', ^/ on the surface of 
the ellipsoid is expressed in ellipsoidal harmonics by 



p 



3M 



+ ^ 



Wi (hD c/ (f ) i/ ii^) c/ {ct>y 



[WM(s:/m'pd 



(T 



By the formula (51) of ^' Harmonics," the external potential at the point v, p, ^ of 
the inequality is 



k 



kp 



lW(i^}^iWJp<^^ 
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But if R is the distance between the point v^^ [x^ <^' on the ellipsoid and the 
external point v, /x, c^, this potential is 1/R, 

If we imagine a particle of mass M/X situated at v^ jjl^ (f>^ the above expression 
multiplied by M/k is the potential of the particle at the point Vq, pi, <^' on the 
ellipsoid. 

We have no need for the general expression for the potential of a particle situated 
anywhere in space at the surface of the ellipsoid, because it is only necessary to 
consider the case where the particle lies on the prolongation of the longest axis of the 
ellipsoid. In this case 

where r is the distance of the particle from the centre of the ellipsoid. 

But it is now no longer necessary to retain the accents to ja', <f)\ since they are 
only the co-ordinates of a point on the ellipsoid. 

Thus the potential of MjX, lying on the longest axis of the ellipsoid at a distance r 
from the centre, at the point vq, [jl, (jy on the ellipsoid, is 

M fr\ 3if^ ®/ft)W(l)®/(i^)f/(^o)W(^)e/(</>) 



,,lPo(.o)©o^y + ^^ 



\3i ii^) (S^i {4>)J V d 



or 



For the types of functions denoted in ^^ Harmonics" EES, OOS, OES, EOS, 
f / (1) = 0, and for EOC, 000, C/ (I-tt) = 0. The only types for which %^ (1) C/ (Itt) 
does not vanish are EEC, OEG; that is to say, cosine-functions of even rank. 
Accordingly the functions left are P/C/ for i and s even, and ^/C/ for i odd and s 
even ; we may however continue to allow C/ to stand for both types. 

For brevity write 

Thence, since ^o (^o) = 1> the potential may be written 



M 



'r\ , s 2*+l -, Jt 



^\lj "^? ^ ®^\l) ^' ^^^ '^' ^^^ ^' ^'''^ '^' ^^^ '^^ ^"^^ 



> for all even values of s. 



It must be observed that p/ and C/ occur as squares in %i ; they also occur 
twice in the numerator in the forms p/ (1) p/ (/x) and C/ (|-'/r) C/ (</>). 

Again @/ ( r ) i^ ^f dimensions — 1 in ||i/, and therefore 0/ ( ^ ) ?P/ (^o) is of zero 

dimensions. From these considerations it follows that p/(/>t) and C/(^) may be 
multiplied by any factors without changing the result, and further that ^/ {y) may 
differ in its mode of definition from ^/ (/x) without producing any change. 

2 A 2 
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The higher harmonics will be considered later, and for the present it is only 
necessary to consider the terms defined hji= 1, s = and i = 2, s = and 2. 

The following are the definitions of the several functions, in accordance with 
" Integrals" : — 

(v) = z., f 1 (im) = /x, Ci {<!>) = v/(l -k'' cos^ </.), 



where g/ = |^[1 + k^+D] and D^ = 1— kV^, with upper sign f or 5 = 0, and lower 

for 5 = 2, 

Hence 

1 (1) Ci (irr) - 1, f / (1) C/ (iTr) = (k^-^/) q/' {s = 0, 2). 



Then from ^^ Integrals," equations (5) and (6), 

2^ 



El = 1, K. 



2 — -^ 
O 



[D'± (1 +iKV^) (1 --2k^^) D] {s = 0, 2). 



Thus as far as the second order of harmonics the potential of Jf/X at Vq, [jl, (f> is 



M 



IK 



kk 



oi^)+3€,(|)fi(^o)fi(/^)CiW + J^ 



1 



+ 



€2' 



{^'-<li)% 



2^'€/(T)f/(^0)i/(/x)C/(<^) 



1 



We must now express the several solid harmonics involved in this expression in 
terms of x, y^ z co-ordinates of a point on the surface of the ellipsoid. 
We have 



z 



f 1 (^0) f 1 (ft) Ci (^) = vi, y/{i--^ cos^ <^) - I • 

By the definition of ellipsoidal co-ordinates the three values of w^ which satisfy the 
equation 



X 



2 1 / 2 ' 2 

0) •— 1//C Ct) 



-j- __-^ — + ^ _ p = are v\ jjl^ -^ (1 — k^^ cos^ <^). 

Ct) — 1 CO K 



Hence we have the following identity 






0) 



2 



1 / 2 ^ 

I/k Ct) 



S^ + ^ __ y^2 ^ P Mzi^) (^'--^') ( 1-^'' cos' (^-coV) 



2 ' 2 i ' 2 

1 0) 



(a)^-l/K')(ct)'-l)coV 



9 
Putting o)^ = i| (s = 0, 2) we find 



L JU 



.2 



1 



2/ 



,2 



i z^ 






'2 P 



K^ — i 



g/ ^^ ^i/ ^ 



2 2 



{s = 0, 2). 



FIGURE AND STABILITY OF A LIQUID SATELLITE. 



181 



Hence the potential of MfX is 



a^«a)-^*©-i 






X 



siBA^ 



€, 



q,%"{K'-qoy- 



5(J|; 



•'i) 



ST/ 



?2V {•<'-q2'y 












E2 
5®2 (r^ 



qo'qo'Hx'-^o') 



'€/ft) 



E2' 



^Uiq^'{>^-qi) 



5®j 



r 



K 



Es 



5(^2 ( r^ 



^qo'H>c'-qoy ^—^q^'i'^-qij 



%i 



,/r\ 



E: 



^g„V («^-?oT+ -^=F^ ?2V ('^^-^27 



21:2' 



For the object immediately in view we only need the terms involving x, y, z, and 
may therefore drop the first and last terms. 

The expressions for g'/ and for Sa^ in terms of k^ have been given above ; by means 
of these I find that 



go" {>^-qo') _ 9g/ g/^ (K^-qJ) __M 



Wi 



ADk 



2' 



Ei 



^JL/K. 



(note the interchange in the suffixes of the q^s). 



2.V 2 



A common factor -4- ^^l may be taken from all the coefficients of x^, y^, z^, and 



since ^^0^2^ 



= ^K^ this common factor is equal to 15/4Z). Hence the terms in x^^ y^^ %" 
inside { } become 



15 x^ 
4i> F L 



T 



T 



- €2 X q^ ix'-qo') + €/ X q^" i-'-q^') 



.h 






+ 



42) F 



®'©*" 



^ 



+ ®/ (| ) ?/' 



/2 



€2(0|j(«^-<?o^)-€/(j)!v(-^-'?/) 



M I q^ 



On substituting for the several coefficients their values in terms of k', I find that 
the potential of if/, -y be written in the following form :- 
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uh4 









{2k'-1) 



(2-/^) 



2 



2 ^tJ^2 



(©2+€/) 



(€2+ ©2') 



+ 



5 ^ 



(i+K^) 



D 



+ 2(^2+®/) 



• ■ 



• (6). 



It may be observed that this satisfies Laplaoe^s equation, as it should do. 
It remains to obtain approximate expressions for the ©'s. 
The expression for these functions is given by 



<a* (^) = f / {^) 



•00 



dv 



4f/WT(^'-Jr(^'-iAT'* 



We require these when v^ which will be put equal to r/k, is large ; thus we must 
develop in powers of 1/v, 

Now 



(y2_iyi/2(^^2_y^2yii2 



V 



i + _ o o + ^o . -. + ^A a ^ + 



2/cV 



2,2 



3_4,4 



2'k^v 



4. ,6. .6 



2*kV 



Since ^o(i') = 1, we have by integration 



©0 



'v\ 



k 



.k 



■ , l + /c^ k\ 3 + 2«:^+3/c* k" , 5 + 3k^+3k*+5k^ ¥ , 



&i^ r' 



40k' 



ry-n 



112k 



6 



r 



,6 



. (7). 



There is no immediate need for this term, since it has been omitted above, but it will 
occur again hereafter. 

Since ^1(2^) = v^ we have 



3©! 



r 






, ^3(1 + k') F ^3(3 + 2K'+3K')k' 



10k" r 



56k' 



+ . 



r 



» » 



(8). 



Lastly, since '^2(1^) = v^—q^JK^, we have 



[f/(^)] 



-2 



V 



1+^ + ^^^* 



^u' 



4 4 

K V 



4- 



(s = 0, 2), 



so that 



1 



V 



6 



1 + 



\±^ + ?^^4 + fi±M+3£' + (1+'^'')^// .. Ht\l 



2k' 



2 2 

K J V 



8/c^ 






K 



4 / 4 



"1 e ♦ * 



(s = 0,2). 



If we integrate this, multiply it by '^'i{v) and write rjk for v, we find. 



^®>'(i)=*; 



1 + 



5(1 + k") 3g/ 
14/c' 7k' 



F 



rp^ 



+ 5 



2 



3 + 2/<:^+3/c^ 5(7/ c/, 
>2k' i26/c'^' 21fc^ 



/c^ 



r' 



.4* 



(6' -0,2) (9), 



FIGURE AND STABILITY OF A LIQUID SATELLITE. 



183 



On substituting for q, its value, I find 



5®/ 



Whence 



r 



„,3 



■ 7(1 + K^) + 2Z> F 5{11 + I0/c"+11k^ + 4(1 + k^)X>} ^ ' 
14k^ r^ 168/c* 



r 



5[©2+€/] = 



2/c^ 



s 



1+k' F 5(11 + 10/<'+1Ik^) ^^ 



D 



[(©a-M 






2yb' 



^^ '2ie r"^ 168/c* 



r' 



J_ F 20(1 + F ) F 
7k' r' 168/c* r^ 



Substituting these values in (6) we have for the potential of Mj\ at the surface of 
the ellipsoid, as far as concerning terms involving «, y, z, 



M 

U 



3©i r .;- 



r\ 


z x^ ¥ 


V 


' k 2F • r«- 




y' ¥' 




2F ' r" [ 



1 + 



3(3 + k') F 5(5 + 2/cM-k^) F 



14/c' 7-2 



56k' 



r 



,4 • 



• 3(1 + 3k') F 5(5k'^+2k'+1) F ■ 
14k' r''^ 56k* r'"' 



z' P 

F • ^3 _ 



• 3(1 + k' )F 5( 3 + 2k'+3k*) F 

7 k' r" 56k* r*'" 



(10). 



If the system be rendered statical by the imposition of a rotation potential, we 
must add to the above such a potential, and that of the ellipsoid itself. 

The expression for the internal potential of an ellipsoid v^ is given in (65) of 
'* Harmonics" ; it is 

2yfc lio (^o) F P,' {v,) ¥ ii^ (^o) F ii (^o)r 



I will now introduce an abridged notation which was used in some of my previous 
papers, as follows :- 

Then, since 

PiM''o)=/\/(''»'-^^)' W(^o) = y(^o^-l), iiW = Vo, 

we may, on omitting the term independent of x, y, z, write this potential in the form 



3M 



X 



Ai^+ 



y 



z' 



2k LF (vo'- 1/k') " F(vo'-l) 



^''^F^'"^' 



'0 



(11) 



The rotation with angular velocity co takes place about an axis parallel to x through 
the centre of inertia of the system, which consists of two masses M and M/X distant r 
from one another. Hence the rotation potential is 



2(0 



y +{^- 



r 



l + X. 






a)V 



(1+X)M h} 2{i+xy ' 



(12). 
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The last term, being independent of x, y, z, has no present interest. Then, collecting 
results from (10), (11), and (12), the whole potential, as far as material, is 



2 



M 
k 



M'Hl 



Xw^k 



27,2 



r 



(1 + X)MJ 



z 
h 



'- 1 Ik') 1 3X " r' L 1 4.1^ ~" '? 56k* r* " ' * J J 



F(l'o^-l//C^) 






1 Vq'-I F 

' 3X r» 



3(3k^-H) F 5(5k*+2k'+1) J:* 



0) 



Vc' 



z 



2 



2,. 2 



k\ 



^' 3X r« 



• 3(l + ;c^) F 5(3 + 2k^+3k *)A;* ' 
7k r 



SGk" 



r 



co^A;^ 2 



3if 



(^0^-1) 



~3M''" 



The condition that the figure of equilibrium should be the ellipsoid of reference is 
that this potential when equated to a constant should reproduce the equation to the 
ellipsoid. The coefficient of z must therefore vanish, and the three coefficients written 
inside { } must be equal to one another. These conditions give the angular velocity 
and equations for determining the figure, but as the subject will be reconsidered from 
a different point of view hereafter, I do not pursue the investigation here. 

At present it need only be noted that the coefficient of z vanishes, and that the 
three coefficients are equal to one another. It is clear then that the potential U of 
the system, as rendered statical, may be written 



U= - 






p^^—1/k^¥ 



3X 



/y%^ 



1 + 



3 (3 + k') F 5 (5 + 2/c' + /c') k' 



14/c 



2 



+ 



T 



56k 



4 



..4 



X i 






.2 






y' 



+ 



z 



[F{v,'--ll^) F(ro'-l) l^Wl 



Now gravity g at the surface of the ellipsoid is —dU/dn^ where n is the outward 
normal to the ellipsoid. 
Hence 

px 



9= - 



a Try ^ TT ^ •jry 

l/ , py aU ^ pz da 

-j- ^ — ~-~ — — . -j- - - ■ 



UMV-l/'c') dx F(fo'-~]) dy kW 8^ 



Now 



X' 






y 



2 



+ 



Z 



2 



4. 4' 



p' ¥{v,'-IIkJ F(V-1)' T<^\ 



and in our alternative notation 



^0—2 
K 



1 _ cos^ y 
sm p 



Therefore 



9 



pk 



^ 3X?-^ sin^ yS 



I — li I r 



3(3 + K^) F _^ 5(5 + 2kVk*) F 



14/c' 



/y* 



2 + 



56k' 



rp^ 



4- 
I • • * 



• (13). 



FIGURE AND STABILITY OF A LIQUID SATELLITE. 185 

As already remarked, we shall put M = ^rrpsb^ - — -; also, since the three axes of 

1 + A 

the ellipsoid are k cos y cosec /3, k cos /3 cosec /3, k cosec /3, we have 

F cos /3 cos y _ Xa^ 
sin^yS l + X* 

Hence 

k^ cos^ y _ Sb^ cos y tan /3 
3Xr^ sin^ ^ 3(l+X)r^ ' 

and the coefficient of the series in the expression for g does not become infinite when 
X vanishes ; however, it is perhaps more convenient to leave it in the form as written 
above. 

This expression for gravity is the result required, but it is to be noted that it is 
determined on the hypothesis that the distant body is a particle or sphere instead of 
being an ellipsoid. 

The development ceases with terms of the seventh order, and the harmonic terms 
of third and higher orders have been neglected. Now the harmonic deformation of 
Roche's ellipsoid of the third order of harmonics is of order k^/r^ in inverse powers 
of r. This deformation is treated as surface density., If we were to proceed to closer 
approximation, we should have to take account of the square of the thickness of the 
layer ; such terms would be of order k^/r^. Since, then, we are avowedly neglecting 
terms of this order, it is no use to carry the development higher than terms of the 
seventh order. 

§ 6. Form of the Expression for the Gravitational Lost Energy of the System. 

The system consists of two ellipsoids, say e and E, with their longest axes co-linear, 
and each of them is distorted by deformations expressible by ellipsoidal harmonics of 
orders higher than the second. To the order of approximation to be adopted, these 
deformations may be replaced by layers of surface density, which may be denoted by 
I and L respectively. 

The lost energy of the system may be represented symbolically by 

F=i(e+^)^+|-(^4-L)^ + (e+/)(^+X). 

Let 5, S denote two spheres of masses equal to e and E and concentric with them 
respectively. 

Then the whole may be written 

V =^ ^ee-\-^EE-\-eE-\-^lU^LL-^ {e-^S)l-\-{E-\-s) L 
-\-{e-^s)L-\-{E-^S)l + lL, 

In the term SI I divide S into two parts, namely /Si, which is to contain all the 
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terms in the potential of S at the surface of e excepting terms expressible by 
ellipsoidal harmonics of the second order with respect to the ellipsoid e ; and 82^ which 
is to contain the omitted terms of the second order. Similarly, the term sL is to be 
divided into s^L and s^L, 

Take the centre of e as origin of co-ordinates x, y^ % with the % axis passing through 
the centre of £?, the y axis coincident with the mean axis of e and the x axis coincident 
with the least axis of e. 

Since I is expressible by harmonics higher than the second order, and since i/^+ 
is expressible by harmonics of orders and 2, it follows that the moment of inertia of 
the layer I about the axis is zero. If therefore co is the angular velocity of the system, 
a contribution to the lost energy of the system which may be written symbolically 

It follows therefore that we may write 

Similarly, if the ellipsoid E be referred to a parallel co-ordinate system X, F, Z 
through its centre, and such that 

cc = X, y -=■ Y ^ z = Z^r, 
so that r is the distance between the two origins, we have 

The problem is already so complicated that it will be convenient to omit certain 
small terms in the expression for the lost energy^ which it would be very troublesome 
to evaluate. 

The term {e—s)L represents the mutual energy of the departure from sphericity 
of e with the layer of surface density L on E, This term is clearly very small and 
will be omitted. Similarly (E—S) I will be neglected. It will appear from the 
results below that these terms are at least of the seventh order in powers of 1/n 
A fortiori IL, which is at least of the eie^hth order, will be omitted. 

The whole expression for V will now be divided into several portions. 

Let {eE\ be that portion of eE in which each ellipsoid may be replaced by a 
particle ; it is, in fact, the product of the masses of e and E divided by n 

Let {eE)2 be the rest of eE, 

Let {vv) denote that portion of V in which the larger body E may be replaced by a 

sphere ; then 

{vv) = lee + ^11 +[6+^2+ i^^ {y^ + ^^)] I + dis- 
similarly, let 

{VV) = lEE^lLL^iE-^s^+^cJ" {Y'-vZ')'\L+.%L. 

Then F= (eJS)i + (e^)2 + {t)t^) + (FF)+ neglected terms. 
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For Roche's problem, when the second body is a particle, F reduces to {eE)i-{-{vv), 
but in the modified form of the problem which I am going to solve the whole 
expression is required. 

The evaluation of (e^) is so complicated that I devote a special section to it. 

§ 7. The Mutual Energy of Two Ellipsoids,^ 

The semi-axes of the ellipsoid of mass e are to be denoted a, 6, c, and the 
corresponding notation for the other ellipsoid is E, A, B^ (7. The distance between 
the centres of e and ^ is 7* ; the axes c of e and C of E are in the same straight line, 
while a, A and 6, B are respectively parallel. For brevity I imagine the densities of 
the ellipsoids to be unity. 

If the external potential of e be ?7, and if d^ be an element of volume or of mass 
of E, the lost energy to be evaluated is 

{eE)^^Uda, 

integrated throughout the ellipsoid E, 

Let us suppose provisionally that the co-ordinates of the centres e and E are x, y, z 
and X, T", Z, and let f, 77, ^ be the co-ordinates of the element (iO ; the axes being 
respectively parallel to a, 6, c or ^, Z?, (7 with arbitrary origin. 

If i2^= {^—xY-{-[r)—yY-^{t,--zf, it is well known that the potential ?7of e at the 
point ^, 7], C is given by 

ex Q I r? TS^ P)^ \^ 1 

U - e'< ^ (r72_£_ 4-/)2_?_ 4.^2_£__\ £ 

^-^t(2n + l)(2n-l-3)2n!\ af H HV R' 
Since ^ satisfies Laplace's equation, we may eliminate ^-^ , and observing that 

;r-^ , ^r-Ts are the same as tt^ , tt-^ respectively, we nave 

pj2 ;^2 ;^2 ;^2 ;^2 



af ay a^^ ' 'd^ ' 'av 

It follows therefore that 






tr-,g (-rs p./i a- a-yi 

o{2« + l)(27i + 3)2ri! Wdx' dyVR' 
Since the operator is independent of f, 7;, ^, we have 



^^ o(2n + l)(2n + 3)27i! VK'aa;^ a^i/V ' J i^ 



^ The results of this section were arrived at originally by a longer method. I have to thank one of the 
referees for showing me the following procedure. — February 26, 1906, 
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dD, . 



But -p- is the potential of the ellipsoid E at the centre of the ellipsoid e, and by 
an exactly parallel transformation 



\~^Et 



{-Y3 



Kl 



2n 



1 a^ 3 



2 \n 



{2n+l){2n + 3)2n\ ^ \K^ dx^ dfj p' 

where p' = (x-Xf + {y- Yf + {z-Zf. 

Since our co-ordinate axes have a perfectly arbitrary origin, we may at once piit 
X = 0, Y = 0, Z = r, z = 0, and after eflPecting the several differentiations put .x^ 0, 

It follows that, on putting x — 0^ ^ = 0, after differentiation and writing 



00 

{eE) = eEt 



(-)"3 



k 



2n 



a d' . d'Y 



00 



+ 7^J X 



i-ys 



'^2 

C 



-XHr^.-^.+ 



aM'i 



o{2n+l){2n+S)2n\ W dx^ dyV o {2i+l){2i+3)2i\ ^\K'dx' dyV p' 

X 92 92 1 3^ 3^ 

If we denote the operator -^ ^--^ + -^-^ by d^ and the operator --^ 3— ^ + ^-^ by D' 



K dx^ dy 



K^ dx^ dy' 



we have 



(eE) = eE 



l_±Pd'+^J^/cid*- 



2.5 



2^5.7 



2\3\5.7 



k'd'. 



2.5^ ^2^5.7^ 2*.3^.5.7 



iC«D«... 



1 



= e^ 



1_ _L(W+ ^'D^)+-^s^{k'd^+ ^^1+ 2^^^ied^D^ 



2*. 3^. 5. 7 



{/c«d« + i^«D«) 



1 



2*. 5^7 



(k'Ktd'D' + k'IQd*!)')... 



Jp 



On effecting the several differentiations, and putting x = 0, y = 0, we find 



d 



2I 



l^i, + l^; d''^=^^J^ + %+3]; Dw'^ = ^J-L> + k + i-o + 3]; 



r^W 



r^ W K- 



r^/c^K" k" K^ 



d«- = 






3^5/ 5 



r' \k* ' k* ' k^ 



r 



12 2 1 

. — + ■ ■■ + — ■ 



/c*K^^K^"^K^K^^/c^"^K^"^^'^ 



and the remaining functions may be found by appropriate changes of small and large 
letters. 

If now we again use p to denote the density of the spheroids, and revert to the 
notation employed elsewhere, namely, 



e = ^ 
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we 



find 



(ei?) = {^trpay 



(1+xy 



1 + 









+ 



+ 



2.5r^ 
3 

3 

2\ 5 V^ 

1 

2*.3.7r' 

9 

2\5.7r' 



P(i„+1^ + ^2^1 



K 



K: 



+ 1 



fQ j — __|-_|-3j-|- J^ 



¥Kl 



4 
K K 

\K XV. 



2 

4- _. 4- Q 






H — ji + rrs + 3 



K 



KI: 



^«(4 + ^ + -^,+ 5) + i-(A + |- + | + 5 



6 4 



K 



Fi^^ 



5 1.2.2.1.. 



/v. XV xv /v. Xx, xv 



+ tf^ 



1 2 

+ 4 + 






.k'K^ ' K* ' k'K^ k' K^ 



(14). 



The first term in this expression is that which was called above {eE)i, and the rest 
constitutes {eE)^* 

If the body E were a sphere, the only portions of (14) which wonld remain would 
be the parts of the expression independent of F^ 

With the object of effecting certain differentiations hereafter, it is desirable that 
the formula for {eE) should be expressed in terms of the semi-axes a, &, c and A,B, C. 

In accordance with the notation used elsewhere, we have 



a = 



k cos y 7 ^ cos B 



sin ^ 



sin^ 



c = 



k 



sinyS 



, where sin ^ = /c sin y, 



^ = /^ cos r ^ ^ i^cos B ^ (j^_jL^^ where sin B = K sin r.^ 
sm B sm B sm B 

The result of the translation into this other notation is as follows :- — 



,.) = «.,.,^4..^p.-.._.._...,.,C] 



+ 



+ 



2^5.7r' 
3 



[S{a*+¥) + 8c''~8(^{a'+¥) + 2a'¥+ssimeinA,B,C] 



2^.5V* 






2*.3.7r 



+ 



9 



2*.5.7r^ 



[2 (^ V + B'¥ + (7 V) + (4^ + 5^ + C) {a' +¥+ c') 

, [1 6c« - 5 (a« + 6«) - 24c* {a^ +h^) + l8c' (a* + &*) 

— 3aV (a^+6^) + 12a^6V+ same in ^, i?, C] 

+ WA'{3a'+¥-4(f)-2A'B^{a'+b'+2c') 
+ same with small and large letters inter- 



changed] 



(15). 



* The fact that capital ^ is nearly the same as B roust be pardoned ; it cannot, I think, cause any 
confusion. 
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§ 8. Remaining Terms in the Expression for the Lost Energy, 

If e be the mass of an ellipsoid of semi-axes a, 6, o, the lost energy of its 
concentration is 



where 



i {ee) = f^e% 






du 



In the present case 



{u-\-aY'' {u-^hY k^-^<^'f'' 

4. a ^ 



and 



/ 2 



oo 



dv 



where 



I'D 



^ = 



Thus 



1 ^ 1 
K sin y sin j8 * 



By symmetry 



where 



1 (ee) = 1% (fTrpa^f ~^~-p^ ...... . (16). 

i(^^) = ,%(|7rpa7^^^ (17), 



/» GO 



^ 



(iN 



and 



KJ^,{l^'-lf\n^-llK'f'' 



n 



K sin r sin B ' 



The lost energy {IS-^ is the potential of a particle /S, equal in mass to E placed at 
the centre of i?, with the omission of terms of the second order of harmonics, 
multiplied by the density of the layer I and integrated over the surface of e. This is 
the same as the potential of the layer l^ with the omission of harmonic terms of the 
second order (and there are none such) at the centre of E multiplied by the mass 
of^. 

A typical term in the surface density representing the layer I is, say, 

mi (/^) €i m 

The external potential corresponding to such a term, at the point ^, /x, <^, is by (51) 
of '' Harmonics,'' 

f (^''''^' ITl) •^''^^ (^>^^^ (^") ^'-^ ^^^ ®^^ ^'^)- 

The co-ordinates of the centre of £ are j^ = y,/x = l,<^=|^7r; and the mass of E is 
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Hence the contribution to (ISi) corresponding to this term is 

I (t'^pa^)' (~X)"^//<®/(^V^)W('^)i/(l)C/(i7r). 

Now |P/{l)C/(i7r) vanishes for all harmonics except cosine-harmonics of even rank. 
Therefore 



(IS,) = ^ iiTrpaJ ^-^ t//®/(#)f /(vo)f /(I) j g^, (in) 



€ln 



% even 
i odd 



. . (18), 



the summation being for all values of i greater than 2, and for all even values of s. 

The lost energy {Ls^ is expressible by the similar function for the other ellipsoid, 
but I have not adopted a specific notation for the ellipsoidal harmonics for the 
ellipsoid E^ and therefore cannot write down the result. 

The lost energy {e-\- S2-^\o^^ {y^ -\-z^))l is the potential of the ellipsoid e, together 
with the potential of S in as far as it involves harmonics of the second order, and 
a rotation potential, multiplied by the density of the layer I, and integrated over 
the surface of e. That is to say it is the potential of gravity on the ellipsoid e 
integrated throughout the layer I, which it is not permissible to regard as surface 
density. 

If the thickness of the layer be ^, and if dl, be a slice of that part of the layer 
which is erected normally on an element dcr of the surface e, then pd^dcr is an 
element of mass of the layer. The potential of gravity is — ^^. Hence the lost 
energy is 

I gC dC dor = -|/) J ^f da-. 

Accordingly the lost energy is equal and opposite to the work done in raising the 
layer, considered as surface density, through half its thickness, against gravity. 
We may take as a typical term 



c=pfmiiH')<^H^)> 



and we have shown in (13) that 



9= - 



X 



pk 



fTrpa 



1 + X 



^1^ ¥ cosy 



'^ ^ 3(3 + /^) F ^ 5{5 + 2k'+k') ¥ ' 
14k^ r^ 56k^ 



r 



It should be noted that this expression for gravity takes no account of the change 
in the ellipticity of e which is due to the fact that E is an ellipsoid and not a sphere. 
The error introduced thus is however outside the limits of accuracy which have been 
adopted. 

Accordingly this portion of the lost energy is 



||7rpV^-^(y:.r{A/ + seri6s} 



mHH-wm'pda. 
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Now in 5 5 we defined 



Hence 



«' = ss/(^f»'<^'«<*ff^''- 



{<!+S,+W+2'))(=-|(|^pa»fp-^,2i££{A.'+seri<s8}e/ . (19)^ 



the summation being for all harmonics. 

In determining the lost energy ^{U) we may treat the layer as surface density. 
A typical term in the surface density is />p/-^p/(/x)C/(^), and the surface value of its 
potential is 



|(l'^pa«)Y^//i/(''o)a/(«'o)W(/^)e/(.^). 



Then, since 



a' = f/(vo)(!|/(''o), 



a typical term of ^(U) is 



Thus 



If iinfM') ^ ( fifuf \ m M «<• m'p dT. 



H«) = |»"paT(^.xi^g/K/. . . 



• . (20), 



the summation being made for all harmonics. 

The value of |-(LX) may be written down by symmetry. 



9. Final Expression for the Lost Energy of the System. 



We have 



V = {e.E\ + M + ( FF) + (eE),. 



The several parts are to be collected from (14) or (15), (16), (17), (18), (19), (20), 
and we have 



{eE\ = (lirpa^ ^-^ 



V 



{vv) = iinpnj 



X 



(1+X)' 



l^X^jj + I tfmi' (r/k) ?p/ (v,) ^/ (1 ) e/ {^tt) {i > 2, s even) 



I o^^ C-gm 



k 2i+l 



(^^^Ai^_^-^«'r /. ■3(3 + K^)F 5(5 + 2k^+k^) F\1 



¥ COi 

SXr^ sin^ ^ 



(all 
harm.) 



(^VV)= symmetrical expression with 1/X in place of X. 
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\ 



{<^E).={i.p^r^^^j[-^^ 



mk+i\ + Kt(~ 



+ 1 



K 






2^5.7r» 



*'(M+^) + ^* 



*' — + — - -I- ^ 



/ J 



+ 









I 3 



2\ 5 V^ 






5. + . 
K 



1 



2^3.7rH 



si 5 , 3 , 3 , p-\ , re/ 5 






"I" rZT 4^ ~~7i "f" 5 



9 

2*.5.7r' 



Fi<c*(i+i. + i+4+4+5 



2t^2 



fC I\. xx, fv XV. Juk. 



■2 ^ 



+ h'IQi i + i + i + -„ + 4 + 5 



2. J 



\KV K* KV K^ K 



(21). 



§ 10. Determination of the Forms of the Ellipsoids. 

We have obtained in the last section the expression for F, the lost energy of the 
system. 

The harmonic deformations of the ellipsoids being of orders higher than the second 
do not enter into the moment of inertia to the order of approximation adopted. 
Hence the moment of inertia about the axis of rotation, which passes through the 
centre of inertia of the system, and is parallel to the a and A. axes of the ellipsoids, is 
given hy __ 



I 



f7r/?a^ 



(22). 



If/ denotes any one of the parameters by which the system is defined, the 
condition that the figures shall be in equilibrium is 

The parameters defining the system may be taken as r, the distance between the 

two centres, cosy, cos ^ for the smaller ellipsoid and cos T, cos B for the larger one. 

Besides these we have the coefiicients f/, F/ of the harmonic inequalities of order i 

and rank s on the two ellipsoids. 

For convenience write 

a = cosy, b = cos/8. 

These letters are chosen on account of the association of cosy, cos/3 with the 

semi-axes a, h of the smaller ellipsoid e. It is unnecessary to adopt a corresponding 

notation for the ellipsoid E^ because, when the problem is solved as regards e, it 

affords the solution for E by symmetry. 

Since 

F cos ^ cos y cosec^ fi = Xa7(l + X), 
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we have 



Hence 
a _h cos y 



a \1 + X/ (abf^ 



X \V3a2/3 



X X^^^r^' G 



h 



X \^/'^ 1 



h _k cos /3 _ 
a^asin/S \l + \J i^^^' a~asin^~U+X/ a^^^' a "~ a sia /3 ~ Vl + X/ (aft)^^^* 



Therefore 



^da _ ^du cfft odh _ da , o ^& 

a a tr 6 a 6 



dc 



da_ # 
a T' 



Therefore 



3a— = 2a 

3a 



da 36 3c 00 



3 oi 3 3 
da CO ac 



(23). 



These enable ns to differentiate, with respect to a? &, functions expressed in terms 
of a, h^ c; the parameters r, fi always occur explicitly. 
The equation of condition for the parameter r is 



dr dr 

On differentiating (22) we have 

-,2 3/ 4 3 Xct)V 



dr ^"^^ (1+X)^' 

In order to differentiate V we must take separately its several portions as defined 
in (21). 

Now 



^ (m) = (f^rpa')'-^ ^ 



dr 



(1+xy 



I ^y;^f /(vo)f /(l)C/(i7r) i®/(yj (^■>2, s even) 



dr 



A 









1 + - 



5(3 + K^ ) F ^ 5(5 + 2 ^^ + ^-^) h 



14:K^ r^ 



56k' 



t 



A 



(all harm.) 



3 



_ __ ( FF) = symmetrical expression for larger ellipsoid. 

or 



dr 



{eE\ = {^irps?)' 



3 



(1+X)n2.5r^ 
3 



f(^, + i)+^^(^3 + i) 



—■■' " I "" " 












2l7r« 
3 

1 

2*.3r^ 

2*. 5r« 



^ 



~4 ' 2 
K K 



+ -. + 3)+^*(| + | + 3 



Fi^M -|r-„ + -„ + 4 + 3 



2 



i„ 



xk'K' k' K 



K K K 



3 






/ mJk 



,,rrAl 5 12 



2 



k'K' ■ K* ' k'K^'^K''^ ■■''^^ 



1 

K' 



+ ¥Kt 



i zJ Zj X 

I ■jiiii->"iiitjiiM fc. L—M ,_ ,. . .^ fcg 1 ■ j^j ^,1 II II I I _ ■ 4>jW mil _ 11^ ||||^,_,ly|^ ^^m^^^mmm 

.KV K^ k'K' k^ K^ 
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orCt) 



The sum of these last four expressions is equal to — 3F/3r, and therefore equal to 
9//8n 



Now let 



o)V = |,r^a^(l + {), 



so that C represents the correction to Kepler's law of periodic times on account of 
the ellipticity of the two bodies. Then we have 



^ = 



lOr' 



y{b^')Mk 



+ 1 



+ 



+ iFi^M i + -4 + -„ + 3 ) + 



56r* 
1 



k 



2 



"l+l+=')+^*(l'+r>+^ 






2 Or* 

9 

8 Or' 



\k'K' K 



2 ' 2 
K 



48r* 



¥ 



3 



K 



,+;^.+-.+5 +i-|^+^,+^+5 



6 



F^«(i+ 1 



+ i + l + i. + sUw' ' 



2. .2 



k'K* K* K'k' K 



K 



, i , 'CI ^ 1 p- 

I 7 1 Ti — n I n I TZli I «^ 



KV K* KV K-* K 



2. .2 



I tMf (-o) f / (1) C/ (Itt) r^ a €/ ft 






14fc' 



r 



(^>2, s even) 
(all harmonics) 



. . (24) 



When ®/ {rjk) is developed in powers of 1/r, its first term is one in r ^'"^^^ ; hence 

r^ ^/ begins with r"*. Now £•* will be determined from terms in the potential of 

or 

the ellipsoid JS of the i*^ order of harmonics, and will therefore involve 7-"^*"^^^ 
Therefore in the series contained in the last term but one of I each term is of order 
^-(2i+i)^ Since the lowest value of i is 3, the term of lowest order in this series is one 
in r~\ and as I shall not attempt to evaluate C beyond r~^ the whole of this series is 
negligible. 

Again, since {/if is of order r~^'~^, and since r~^ occurs as a factor, each term is of 
order r~^*"*^ Thus the lowest term is of order r~^^ and is negligible/^ 

It follows that the only sensible part of C arises from the portion of V denoted 
{eE)2, and the last two terms of (24) may be erased. 

We next consider the parameter //, and, since / does not involve it, the equation 
reduces to 9 V/dfi = ; or, since V only contains // in the part denoted (w), it 
becomes 9 {vv)/dfi = 0. 

This gives, for ^>2, s even 



y? = -i(2^+i) 



XE/ 



^A^^ 



h^ cos^ y 



3Xr^ sin^ fi ^ 



■ 3(3 + K^)F 



,.2 



(25). 



Since this formula contains X in the denominator, it would appear at first sight as if 

* It is proper to remark that the terms retained in ^ are really of higher orders than they appear to be, 
I recur to the neglected portions of ( hereafter in § 23. 

^ C ^ 
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fi became infinite when X = 0. But this is not so, because when (%i{r/k) is developed 
the first term of the series is one in {k/rf'^^ ; now F = - — -^ a^ sin^ /3 sec j8 sec y, and 

1 "T" A 

therefore the formula for ff involves the factor ( - — -• ) . - or ( ) . . 

^ \1 + X/ X \1 + X/ 1 + X 

We see then that // vanishes both when X = and X = oo. 

This factor is a maximum when X = ^{i—2). Therefore we should expect, ecBteris 
paribus, the third harmonics to be most important when X = |^, the fourth, when 
X = f , the fifth when X = 1, and the higher harmonics when X is greater than unity. 
This prevision is partially fulfilled by the numerical results given below, but it was 
not to be expected that it should be exactly so, because the other conditions, are not 
exactly the same in the solutions for various values of X. 

The formula shows, as stated above, that fi'm of order r'"*""^ The series in the 
denominator affects the result but slightly and might be omitted, except, perhaps, in 
the case of the third zonal harmonic. For all harmonics other than cosine -harmonics 
of even rank j^-* is zero. 

It is now possible to eliminate fi from {vv) by substituting for it its value. These 
terms in (/ov) become, in this way, equal to 



X i 2t' + l 



n|)f/(^o)i/(i)c/(M 



(t ^/>aT ^-^_^^y 1^ 2a %i [^/_- A/-->series] 



2-\ 



> . . (26). 



-J 



When i = 3, this term is of order r""^, and is negligible ; hence we need no longer 
pay any attention to the inequalities on the ellipsoid. However, the formula (25) is 
important as rendering it possible to evaluate the inequalities. 

Since for all inequalities, excepting cosine-harmonics of even rank,/' only occurs 
in the energy function as a square, it is in these cases a principal co-ordinate, and 
{%i — hi — series) is a coefiicient of stability. 

But the like is not true for the cosine-harmonics of even rank, because, when we 
consider, for example, the harmonics of the third order, we see that dW/df'dr is of 
the fifth order smd dW/df'dn, dW/dfid'b {s — 0, 2) are of the fourth order. 

It is clear that the inequalities on the ellipsoid E are determinable by symmetrical 
formulae. 

We must now turn to the equations of equilibrium for the parameters a and b. 
Since differentiation with respect to these parameters is effected most conveniently 
by means of the formulae (23), the portion of F called {eE)2 should be written in the 
form (15). After effecting the differentiations it is, however, best to revert to the 
notation involving k, k, y, i^, K, T ; but as an exception to the general rule as to 
notation, it is most convenient to retain the differentials of F(l + l//c^) and of 6^+c^ in 
the forms involving a, 6, c. As the algebraic processes involved are rather long, I 
simply give the results^ as follows :^ — 
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On picking out the numerical coefficients of the several terms in {eE)2 as given in 
(14) or (15), we see that 



3a ~ {eE), 



_ (|7rpayX Jl(ii), 3(iii) , 3(v) , 1 (iv) , 9 (vi) 



+ 






5 + 



36^(.^).= 



(1+X)' SXhT^ 2\5.7r'' 2l5V^ 2*.3.7r' 2*,5.7r^ 



. . (ii)' . . (iii)' . . (v)' . . (iv)' . . (vi)'. 



• • 



Observe that -^r-r 



2 



¥ 



K^ sin^ y sin^ ^ 



= c^, and write 






3 K^Bixry 



(iii) + 



K^sln^y ^.^^^ _^ 9 K^sin^y 



-(1^) + 



(vi) 



> (27). 



<r 



= . . . . (iliy ..... (v)' (iv)' (viy 



Then we have 



3a^(e^). = (|^paT ^ 



(1+xy 

X 



^ (5^-2a'-2c^)+-i^.Tc' 



86 1 (e£). = (|,paT ^j ^ 



5^3 



5r 



3 



% 



V- 






5r' 



. (28). 






The terms in V denoted {eE)i and (^VV) do not contain a, &, c, and their 
differentials with respect to % b are zero ; also, after omission of the terms in ff, 
(vv) is reduced to 



Hence 



Now 



(1+X)' 



m ^ ^ ■ (l + X)^L \ 3a 86 3c /^ 



(1 + X) 



86 8a 8c 



«I=-|p^^wq^^w=-|a>' 



86 



8t|r _ 

3c 






|fi(»'o)ai(vo) = -|^>. 



Since i/f is homogeneous of degree — 1 in a, b, c, the sum of these three is equal to 
-1, so that 
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Therefore 



3a|;W=-(l^/-aTj^.| 



(3A,'-PW. 



3*^H=-(f.paTpA_,|(3a.'-i^), 



> 



J 



* » 



(29). 



On adding together (28) and (29), we find 



dV 



3a-^ = (|^paT^^:^L 



I (3 Ai^-P«A)+ ^3 (&^-2a^-2o^+ro^)' 



3b 



ay 

d\ 



= (i"-""')"?!^ 



-|(Sa.'-»)+ ^{a"-2(.'-2c'+,Tc')' 



By means of (i) and (i)' we find the differentials of the moment of inertia /; 



ley are 



3aS=-|^/oa«-^.|(6^+c^), 



c?a 



i+x. 



rlT 

3i,|= |.,a3 



^ .|(26^-c^). 



1+X 



Then, since |-<u^ = ^Trpa'. —3 (1 + 0> 



3,. 2, 



,.-a|=-(|.,aTp^.-^(l+X)(l + £)(5-+o-), 



fw% 



dl 

dh 



= il^P^J (~^. • ^3 (J + ^) (1 + (2&^-c^). 



Now the equations for equilibrium for the parameters a and b are 



Therefore 



rta aa 



;!; 






■% 



3X(3Ai'-P»|')+ -3[6^-2a^-2c^+TC^- (1 + X) (1 + {b' + c')] = 
-3X(3aii-p.|»)+ -3[a^-26^+2c^+o-c^+(l+X)(l + 0(26^-c^)] = 



> 



(30). 



y 



Subtracting the second of these from the first and dividing by 9X, we have 



k 



^i'-Ai^ = 5V3K-^'+^Ml+^)(l + 0-*«MT-a-)]. . . . (31). 



3Xr 
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Since we may write 3 Ai' — ^^ct// in the form — ('^i' — Ai') — {^i — Ai'), the first of 
(30) in combination with (31) gives 



^i-A/ 



k 



SXr 



[a' + 2c'+c'{l + \){l + C)-lc'{2T+(r)] . . . (32). 



Referring to the values of r and a in (27), I find 
--HT-c7) = A5^'[2-3(l + K^)sin^y] + /e-^'[4(l + K^)-(5 + 6K^+5K0sin^y] 

■ IK. I K. 






+ il) 717-2 [2K" cos' j-k'^ sin' y ( 1 + 3K')] 



+ 



r K 



9 



[4K" (1 +K') cos' y-K'^ sin' y (1 + 2K'+ 5K*)] 



+ ^Og- .|!^ [4(1- k'K') cos' y - k" sin' y ( 1 + k' + K' + 5/c'K')]. 



rVK 



-H2T + 0-) = ^3^ 4^[7 + 5K'-(3 + K')sin'y] + /o:|^,[ll + 6K'+7K*-(5 + 2K' + /c*)sin'y 



'^^ 0^2 2 



4 4 

r K 



+ -u.-E[7 + 5K'_(3 + K')sin'y] 



+ 



5 ,,.•4^4 



[ll + ()K' + 7K'-(5 + 2K'+K") sin' y 



+ -Sg-^2[11 + 3K' + 3K'+7K'K'- (5 + K'+K' + /c'K')sin'y] 

/ K. J\. 



. . (33). 



In all the cases which we shall have to consider the first of these expressions is 

small compared with the second, because k is nearly equal to unity and k' small, and 

because cos' y is also rather small. 

Now let 

e= _x(^_o.) + ^(|+x)cos';8l 

V (33 Im). 

■,) = -■i(2T+cr) + n'+^) J 

Then, since a = c cos y, h = c cos /3, the equations (31) and (32) become 



kc' 



^i' - Ai^ = .-3 (cos' y + X cos' ^ + e) 



3\? 



a,-Ai' = |^,(3 + )i+cos'y + r/) 



> 



(34). 



J 



Eliminating kc^/'S\r^, we have 



I,-Ai^)(cos^y + Xcos^i8 + e) - (^i^- Ai^) (3 + X+ cos^ + t^) . . (35). 
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This is the equation to be satisfied by the axes of the ellipsoid. If we treat e and 77 
as zero, it is the same as that found by Roche. ^ 

^ The form of this equation is so unlike Roche's, that it may be worth while to prove the identity of 
the two. 

Roche writes his equations in the form 



st (t - s) 



udu 



_ ^(1"^) 



^ t{l-t) 

{3 + X)t-Xs}o{l + su) (1 + kt) R s + 3 + Ajo(l+5t^)(l+ti)E t-\-X 



i 



udu 



05 



udu 



0(1 ■\-u)(l-\-tu),B 



where E^ = (1 + it) (1 + su) (1 + tu\ and s is the square of the ratio of the least to the greatest axis, and t 
the square of the ratio of the least to the mean axis. 
In my notation 



a^ o . aP- cos^ y 



If we write us+l ^ 5HL-2, and change the independent variable from u to ;/', we find 



Roche's first integral 



2 cos^ 13 
sin^ y cos^ y 



V sin^ \f/ (sin^ y - sin^ \f/) 

w 



dip 



j> 



jj 



second 



third 



)j 



»j 



_ 2 cos /5 {y sin^ ^ (sin^ y - sin^ \{/) y, 
sin^ y cos y Jo cos^ i/'A 

__ 2 cos^ /? fv sin^ ;/' (sin^ y -- sin^ ^) j , 



where A^ == 1 - k^ sin^ xj/. 
The coefficients are 



sm*" y cos y J 



cos^ \^ A^ 



• • 



. . (A), 



st{t~s) 

{3 + X)t-Xs cos2 I3[{3 + X)-X cos^ /5] 



cos^ y sin^ /5 



g(l-g) 
S + 3+A 



sin^ y cos^ y 

34-^ + C0S2 y 



k'2 sin^ y cos^ y 



h 



• • • • 



• « 



. . (B). 



Hl-i) ^ . 

t^X cos'^ ^ (cos^ y + A cos^ /^) 

Then Roche's equations are equivalent to 

1st of (A) X 1st of (B) - 2nd of (A) x 2nd of (B) = 3rd of (A) x 3rd of (B). 

But the two equations are not independent, and I will only pursue the consideration of the form 
involving the 2nd and 3rd of (A) and (B). 

Now 

sin^ xb (sin2 y _ 31^2 ^) ^ ^-^^ ^ tan^ ip 

'I I " ' ' — ' ' ' , — L»Oo y » 

cos^i/^A A .A 

sin^ ijy (sin^ y - sin^ xp) ^ cos^ /^ sin^ jy _ cos^y tan^ 1/^ 

C0S2 1// A3 k'2 ■ A3 " /c'2 A ' 

and I have proved in (25) of the '' Pear-shaped figure, &c." that 

K [y sin^ if/ 



x^i 



li®] 



sin^yjo A 



#5 



Ai = 11 VI r— 5 I r — 

^ sm^ y Jo A 



(ii/', 



sm^ V a*^ 



sm^ y 



TOL. COVL — ^A, 



AM ■^•^ 
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It is possible to express this equation in terms of elliptic integrals and to use 
Legendre's tables for finding the solution, but the method is very tedious, and after 
finding a few solutions in that way I abandoned it. It may, however, be worth 
while to mention that 



i^l r\i 



K 



^/-Ai^ 



sm y 



K 



sm y 



1 



K 



^ sm y cos y cos 



2 



K 



-^ sm y cos y cos 






rs _ COS^ /3 p 2 COS^/3-K^^ 

P 2 -^ ' 2/2 

AC K K 



E 



whe 



re 



jF = 



d\ls 



0^(1— K^sin^^//) ' 



^ = 



v/(l--/c^sin^i//)(ii//. 



When the forms of the ellipsoids have been determined, the radius vector becomes 
determinable from either of the equations (34). 

The conditions that the internal potential of an ellipsoid satisfies Poisson's equation 
and that t/f is homogeneous in a, 6, c of degree —1, give the two following equations : — 



a 



h 



^ c^ ahc 



Aji + ^/ + ^1 - ^^ = 0. 



Our two equations for 1/r^ may be written 



k 



Ai^-a/ +-^K+6^\+c^€) = 0, 



3\r 



Aji -% + ^^[{Z+\)G'+a'+c'7j] = 0. 



These four equations afford a determinant by which Ai\ ^i\ ^i may be eliminated. 
On reduction we find 

j^ . Sfaho 

1 ^"^ l/a'+l/h'+lfc' 



Therefore Roche's second integral is equal to — r-s — - — (^i - Ai^), and his third integral is equal to 



2 cos^ p 



K sm^ y cos y 



KK^ sin^ y cos y 



(^1^ - All). 



Using these transformations of the second and third of (A), and dropping redundant factors, we get 

(cos2 y + X cos2 p) (Si -. Ai^) = (3 + X + cos2 y) (^ii - Ai^). 



This agrees with the result in the text when c and t] are neglected. 
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On putting a ^ c cos y, & = c cos /S, and noting that i/; — 



2F 



-—— G^ cos y cos ^ = a^, I find that 



G sm y 



and that 



a 



3 



3. 

2 



6 



2F cot y cos ^ - 2 JTS 

1 + sec^ y + sec^ /3_ 



6 



where 



r^ 3 — cos^ y + X ( 1 + cos^ /3) 

1+X 1 + sec^ y + sec^ j8 



+ S 



. (36),^ 



S = 



1+X 



J 1 + sec^ y + sec^ ^ J * 



However, this is not practically the most convenient form from which to compute 
the distance between the two ellipsoids. 

■^ It is by no means obvious how this formula is consistent with results which we know by other means 
to be true. In the case when X = oo we have a liquid planet rotating with the same angular velocity 
as an infinitely small satellite revolving in a circular orbit in its equator. 

Let us first consider the value of (. In the present case the semi-axes A, B, G pertain to the infinitely 
small satellite, and are therefore negligible compared with terms in a^ b, c. Since the axis denoted by c is 
that coincident with the satellite's radius vector, and since the equatorial plane of the planet must have a 
circular section, we have c = b. 

But since ^ = c cos /?, it follows that /i? = or k sin y = 0. Now y does not vanish for a = c cos y, and 
a is the polar semi-radius of the planet ; therefore k = 0. 

If we consider the formula (24) for ^, expressing, however, the several terms in the form of (15), we see 
that for X = oo 



C =. -A.{G^^a^) + ^, . 3 {G^-a^Y+ ~l^, . 5 {c^-a'Y ... , 



whence 



10r2 



3c2 



56r4 



48^6 



9c4 



5c6 



^ = TT^ ^^'"'y+Kozi «i^' r + tqTg «i^' r 



10r2 



56? 



48^6 



• . . « . i 



i , i i i t, , V /* 



The factor of correction to Kepler's law of periodic times for a small satellite revolving about an 
oblate planet, whose equatorial radius is c and whose eccentricity of figure is sin y, is 1 + (^ where ( is 
expressed by the above series (a). 

Now considering the formula (33 bis), we see that for X = oo and /3 =^ 



1 + A 



Ccos2/^=t> 



V 



1 +x 



C^ 



In (36) we therefore have 



8^- 



1-hA 



^"f € 



3(r) + e sec^ 13) 
1 + sec^ /3 + sec^ y_ 



-2C 



I- 



__ 2( tan^ y 
3 -f- tan^ y J 3 + tan^ y 



When /c = 0, the elliptic integral F is equal to y ; thus (36) becomes 



a^ 



2 



2y cot y - 



6 



3 + tan^ y_ 



2- 



6 



21 tan^ y 
3 + tan^ y ' 3 + tari^ y 



+ 



2 D 2 
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§11. Solution of the Equations. 

In all the ellipsoids of which we shall have to find the axes, it happens that K^^tan^y 
is fairly small compared with miity. Hence it is possible to expand A in powers of 
that quantity. 



We have 



and 



A^ = 1 — K^ sin^ y — ^^^^ 7(1 + kI^ tan^ y). 



■i = — i- \\-U^ tan^ y+M'^'' tan^ y--X 
A cosy '2.4 

-i. = — ~ [1 -\^^ tan^ y + ir^^ tan^ y- . . .]. 
A^ 0,0^^ y^ ^ ^2.4 /J 



From this we easily obtain 



(0^ 0/1. ^\ a^ 



I (1 + 1) fL. = cot^ y [7 (3 + tan2 y) - 3 tan y]. 



27rp "" r 

This is the well-known formula for the angular velocity of Maclaurin's ellipsoid. 

It should be remarked that (35) is identically satisfied by X = oo, k = 0, for when we use the above 
values of € and ?;, the equation becomes divisible by 1 + 1« 

Since { is a symmetrical function of a, &, c and A^ B^ (7, it follows that { is the same in form for X and 
for 1/X. Therefore when we consider the case of X = 0, the formula {cC) gives the required result, but 
c and y refer to the large body which is throughout most of this paper indicated by capital letters. 

Thus for X - 0, 

i = ^ Bin^r+VV (5 Bin^ rJ + ^(^ sin^ rj... 
3 a^ sin2 r g a^ / sin^ V y g a^ / sin^ V \3 

TO y.2 eosW3 p + 5 6 ^,4 (^^^^2/3 p j + T8 ^.6 (^^Qg2/8 p j . • • • 

In the case of X = 0, ^ vanishes ; K^ also vanishes and so also does the angle B. Hence we have 

e = f cos2 /5, ^ = C. 
With these values, equation (35) becomes 

(ai"Aii)[cos2y+(X + 0cos2j8] = (aii-Aii)[3+(X + + cos2y]. 

Hence f plays the part of an augmentation to X. 
With X == the equation assumes the form . 

(i(i-Aii)(cos2y + Ccos2^)-(aii-AiO(3+i+cos2y) . . ..... (5). 

It follows therefore that an infinitesimal satellite revolving about an oblate planet, whose rotation is the 
same as the revolution of the satellite, is very nearly identical in form with a small but finite satellite 
whose mass is a fraction of a s;pherkal planet expressed by {". This curious conclusion follows from the fact 
that if we take equation (35) and put e and 7] zero (which corresponds to a spherical planet and small 
sateUite), we get exactly the equation (b) just found, only with X in place of {. 
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Now from (25) of the " Pear-shaped figure," 



siiryJo A • sm^'y Jo A sm'' y Jo A^ / 

When the A's under the integral signs are expanded, all the terms of the series 
involve integrals of one of two types. If we write 



a = 



the types are 

p sin'> 7 
^ ( (iy = siny 



sm y 



log 



1 + sin y 
COS y 



1 



2n-2 



tan 



2n-2 



y- 



2n-l 



(2n-2)(2n-4) 



tan'^''-*y-|- 



+ ^ ^ (2n-2)(2n-4)...2*'''' ^+^ ^ (2n 



^)---^ ro_ 



2). ..2 



(fl-1) 



f 



^2n 



sm y 7 

ocos^^'^^y ^ 



1 



tan''^ y— 



tan'^^^-'y^- 



2n-l 



tan 



/ Y (^^^"^1/ (2n — 3)...5 ^^^2 / w+i (271— 1).. .3 /£ii _]^\ 



27i(2n-2)...2 



2n(2?z-2)...2 



As it is not quite obvious what interpretation is to be put on these formulae for the 
smaller values of n, I may mention that when ^^ = 0, 1, 2 respectively, the first 
integral is sin y ; sin y (ll— 1) ; sin y [|- tan^ y— f (O — 1)], and the second is sin y II ; 
sin y [^ tan^ y - 1- (O - 1 )] ; sin y [^ tan^ 7 "" iTa *^^^ 7 + 172 (^ "" ■'^ ) J' ^^^ larger values 
of n the interpretation is obvious. 

If we use these integrals and write 



^ m:-A:) = k" [o-o-cry^+o-y*-o-3/c'«...] 



K 



sm 



K 



1 (^i-AiO = T,-T,K" + T,K"-r,K'K.. 



Y 



(37), 



we find 

o-o = i [2 tan' 7 + 3 - (3 + sin' y) O] 

o-i = 3% [I tan* y-f tan' y- 5 + (5 + sin^ y) O] 

0-2 = t^I^CtI tan«y-f tan*y+-3-tan^r+7-(5' + sin'y) O] 

0-3 = 4W6 [fftanV-Ttf tan«y+!f tanV-4tan'y-9 + (9 + sm'y)«] 

To = i[-3 + (3-sin'y)0] 

1-1 = -A [I tan' y + 5 -(5 -sin^ y) XI] 

T. = 1^2% [t5 tan* y-l tan' y- 7 + (7 -sin' y) 12] 

T3 = -.iol^hJ-o^tan^y— A-tanV + 2tan'y+9-(9-sin'y)ft] 

T4 = ^4^ y^ tan« y — i%^5 tan« y + 1 tan* y- f tan' y- 11 + ( 11 -sin' y)0] . (38), 
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It would not be difficult to find the general expressions for these functions, but it 

does not seem worth while to do so. 

The equation (35) for determining the form of the ellipsoid involves the factor 

cos^ y + K cos^ ^8 + e ; if we write 

jf^ Xsin^y 

{H-X)cos^y + e' 

this factor may be written in the form [(1 + X) cos^y + e] [l + if/c'^]. Hence the 
equation (35) may be written 

[ro-Tifc'^+ry^-Tsfc'^..] [(l + X) cos^ y+e] [1+1/k^^] 

= fc'^[a•o-o-lK^^ + o■2K:^^-(r3K:^^..] [3 + X + cos^y+7;]. 
If now we put 



we have 



t^O J Vi — —J V2 Vi, V3 — — V2~ Vi 

o-Q ctq To o-q To cTq (Tq Tq (Tq ctq 



> 



Tn TiK + T2'C ... n I /2 /4 I /6 1 

cr^)— ctik: +cr2K: ... 



Hence our equation may be written 

Vo[(l+X) COsV + e] (i+^y.^^^,^.^.,ye^^^) (1^^^.) ^ ^. 

3 + X + cos^y + 77 ^ ^ /\ / 



Whence on writing 



r --. S + X+cos^y+Tj 
~ (l + X) cos^y + e 



K. — ' — — — — — — — ' 



The determination of L for given value of y involves that of 77 and e, and these can 
only be found from an approximate preliminary solution of the whole problem. But 
when L is known approximately, the solution of (39) is very simple, for we first 
neglect the tei'ms in k!^ and k!"^ on the right-hand side, and so determine a first 
approximation to k'. As a fact I have not included the term in k!^ in my computa- 
tions, because it would not make so much as 1' difierence in the value of cos"-^ k'. 

For Roche's problem when e and 77 are neglected the solution is very short, but 
when these terms are included the computation is laborious. 

We now turn to the determination of the radius vector. 

We have 

ci' ^ 2 ^^ J F cosy cos 5 Xa^ . 

bmce c^ = -on ^^d . % ^ ^ := r- , we have 

sm^ /3 sm^ yS l + X 

kc^ __ 1 Ksiny a^ 

3Xr^ ~ 3 (1 + X) cos /3 cos y ' 7'^ * 
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Therefore 



Whence 



Biny^ ^ 3(l + X)cos^cosy r^^ ^ '^ 

a^ _ Q /I , \ \ cos /S cos y To— Ti/c'^+rg/c'^— Tg/c'^. . . 



sm"y 



S + X+cos^y+r? 



(40). 



Thus a table of values of tq, t^, rg, T3 enables us to compute r for a given value of y, 
when K^ has been found. 

In the following tables the vh and t s were computed for the even degrees of y and 
interpolated for the odd degrees. These functions are found as the differences 
between large numbers, and therefore great care would be required to determine them 

» 

Table of Auxiliary Functions. 



r- 


Log Vq. 


Vi. 


Log vi. 


Log V2, 


Log V3. 


To. 


Log Ti. 


1 

Log T2. 


Log T3. 




30 


9- 


94493 


0^14279 


9^15470 


7^9576 


7^15 


0^010592 


6-85440 


5^9805 


5^222 


31 


9- 


94085 


0^ 


15461 


9-18925 


8^ 


0264 


7' 


25 


0^ 


012124 


6- 


94587 


6- 


1063 


5- 


378 


32 


9- 


93661 


0^ 


16722 


9-22329 


8^ 


0941 


7- 


36 


0- 


013824 


7' 


03507 


6' 


2288 


5- 


533 


33 


9- 


93222 


0^ 


18065 


9^ 


25684 


8- 


1607 


7- 


47 


0- 


015702 


7' 


12194 


6' 


3478 


5- 


685 


34 


9- 


92764 


0^ 


19493 


9^ 


28988 


8- 


2263 


7- 


58 


0' 


017774 


7- 


20666 


6- 


4639 


5- 


834 


35 


9- 


92285 


0- 


21006 


9- 


32234 


8' 


2906 


7- 


67 


0- 


020052 


7' 


28937 


6- 


5780 


5' 


979 


36 


9- 


91786 


0' 


22613 


9^ 


35436 


8' 


3541 


7' 


76 


0' 


022554 


7' 


37022 


6' 


6899 


6' 


122 


37 


9- 


91269 


0' 


24327 


9^ 


38608 


8- 


4173 


7' 


87 


0' 


025294 


7- 


44939 


6- 


7994 


6' 


262 


38 


9' 


90731 


0' 


26150 


9^ 


41747 


8' 


4800 


7' 


97 


0" 


028290 


7' 


52697 


6- 


9069 


6- 


400 


39 


9" 


90167 


0- 


28091 


9^ 


44857 


8" 


5418 


8 


•06 


0" 


031559 


7 


•60306 


7' 


0127 


6" 


536 


40 


9- 


•89578 


0" 


30160 


9^ 


47943 


8- 


6031 


8 


•14 


0" 


035121 


7 


•67777 


7' 


1167 


6- 


669 


41 


9 


•88965 





•32368 


9- 


51011 


8- 


•6646 


8 


•23 


0- 


•038994 


7- 


•75120 


7' 


2192 


6' 


802 


42 


9 


•88327 





•34724 


9' 


54063 


8 


•7259 


8 


•32 





•043203 


7 


•82342 


7' 


3202 


6' 


932 


43 


9 


•87665 


0- 


•37244 


9- 


57105 


8" 


7869 


8 


•42 





•047768 


7 


•89452 


7' 


4200 


7' 


061 


44 


9 


•86974 





•39937 


9' 


60138 


8 


•8476 


8 


•51 





•052713 


7" 


96459 


7' 


5186 


7' 


189 


45 


9 


•86257 





•42823 


9- 


63168 


8 


•9081 


8 


•60 





•058064 


8 


•03371 


7' 


6162 


7' 


316 


46 


9 


•85509 





•45917 


9- 


66197 


8 


•9685 


8 


•69 


0" 


063847 


8- 


•10195 


7- 


7128 


7' 


442 


47 


9 


•84729 





•49239 


9' 


69231 


9 


•0290 


8 


•79 





•070093 


8" 


16939 


7- 


8087 


7' 


567 


48 


9 


•83915 





•52812 


9' 


72274 


9 


•0896 


8 


•87 


0- 


•076830 


8 


•23608 


7' 


9036 


7' 


691 


49 


9 


•83069 





•56656 


9' 


75325 


9 


•1505 


8 


•96 





•084093 


8*30210 


7- 


9981 


7' 


814 


50 


9 


•82187 





• 60804 


9' 


78393 


9 


•2117 


9 


•05 





•091916 


8-36752 


8 


0920 


7' 


938 


5i 


9 


•81263 





•65278 


9' 


81477 


9 


•2733 


9 


•15 





'100336 


8 


•43241 


8 


•1855 


8 


•061 


52 


9 


• 80300 





•70121 


9- 


84585 


9 


•3353 


9 


•24 





•109S94 


8 


•49682 


8" 


2786 


8' 


183 


53 


9 


•79298 





•75372 


9' 


87721 


9 


•3976 


9 


•33 





•119134 


8 


•56082 


8 


•3716 


8 


•306 


54 


9 


•78252 





•81075 


9 


•90889 


9 


•4606 


9 


•43 





•129601 


8 


•62447 


8 


• 4643 


8 


429 


55 


9 


•77159 





•87288 


9 


94095 


9 


•5245 


9 


•52 





•140845 


8 


•68785 


8 


•5570 


8 


•552 


56 


9 


•76016 





•94064 


9 


•97342 


9 


•5893 


9 


•62 





•152919 


8 


•75103 


8 


•6500 


8" 


•676 


57 


9 


•74819 


1 


•01471 




•00634 


9 


•6550 


9 


•72 





•165883 


8 


•81407 


8 


•7431 


8 


•801 


58 


9 


•73565 


1 


•09590 




•03977 


9 


•7216 


9 


•82 





•179801 


8 


•87703 


8 


•8367 


8 


•926 


59 


9 


•72252 


1 


•18517 




•07378 


9 


•7893 


9 


•92 





•194740 


8 


•93999 


8 


•9306 


9 


•052 


60 


9 


•70874 


1 


•28362 




•10844 


9 


•8582 




•02 





•210779 


9 


•00303 


9 


•0252 


9 


•179 


61 


9 


•69426 


1 


•39268 




•14385 


9 


•9287 




•13 





•227997 


9 


•06626 


9 


•1206 


9 


•308 


62 


9 


•67910 


1 


•51349 




•17998 




•0006 




•24 





•246485 


9 


•12966 


9 


•2168 


9 


•439 


63 


9-66319 


r6480 




•21696 


•0743 


•36 


0^266343 


9-19331 


9^3140 


9^572 
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with a very high degree of accuracy. The differences of the tabulated numbers do 
not run with perfect smoothness, showing that there are residual errors of one or two 
units in the last place of decimals. The accuracy is however amply sufficient 
for the end in view, and it would have been wasteful to spend more time over 
the computations. 



§ 12. Determination of the Form of the Second Ellipsoid, 

The parameters r and K determine the form of the second ellipsoid in the same 
way that y and k determine the first. It is obvious that a/r is determinable in two 
ways, and therefore any given value of y must correspond to a certain definite value 
of r. The fitting together of the two solutions can only be effected with accuracy 
by interpolation, but it would be so enormously laborious to find by mere conjecture 
the region in which to begin calculating with assumed values of T, that an approxi- 
mate solution of the problem becomes a practical necessity. 

After various trials I find that on neglecting € and r} and writing 

___ sin^y 

^ ~ 7(3 + X+cosV) 

the solution for k may be written approximately in the form 

^"^ IV'x+rJ-y [1^(9-^7X) x-(lll + 228X-^49X^)x- 

-Tr(^875 + 28185X + 42333X2-3773X^)x'.*0 . . (41). 

If e were added to the numerator of the factor outside the bracket, and r) to the 
denominator, this formula would give nearly as good results as the more accurate 
method of the last section, 

xxiso X unci 

r^ 3 + X + cos^y l -5 \ /a «\ /a 

+ _J__(3533389 + 3222607X+1021479X^+60025X'^ . . (42). 

In order to obtain the desired approximation, it is necessary to express aP/r^ hj a 
series which can be inverted ; but this is not possible in the form just given, because 
cos /3 depends on k and therefore involves x- I find then by means of the above 
series for k'^ that 

cos^cosy= l--|(74-X)x+|{69---106X~7X')x'+i^6{253 + 753X--1901X'-49X')x^./ 



FIGUEE AND STABILITY OF A LIQUID SATELLITE. 209 

On introducing this in the above formula for a^Jr^, I find 

^; = j^(i+x)x[i-(2-x)x-(io9+6m-ox^)x^ 

- 2^^^(157n2 + 261395X + 97656X^+1568X=')x'...]- 
On writing 

^ = 14 



1 + Xr^' 



and inverting the series we find 



^ = a + (2-X)a' + (117 + 59X + 2X')a^ 

+__X_(306872 + 269975X+5739X'+908X^^ . . (43). 

This series expresses a function of y in series proceeding by powers of d^jr^^ and a 
similar series must also connect a function of T with the radius vector, so as to 
determine the figure of the second ellipsoid appropriately. This second series may 
be written down by symmetry. 

Since X must now be replaced by 1/X, the function corresponding to a is ~^^- — ~ — 

• • X sm^ r 

or Xa. and the function corresponding to v is' ^, ^ — r-~ — — = . 

^ B A 7[(3 + cos'r)X+l] 

If then we write 

sm JL 






7[(3 + cos'r)\+l]' 



the symmetrical series for the other ellipsoid is 

X = a + (2X-l) a^ + (ll7X'+59X+2X') a^ 

+ ^,3^^ (306872X''+269975X'+5739X^+908)a*... 

Now a is easily computed for the first ellipsoid, and then X is computed by the 
series. Thus we have 

We obtain in this way a fairly accurate value of r corresponding to the value of y 
which determines the first ellipsoid. We can then compute K^ by the method of the 
last section. We may thus obtain a good idea of the values of r and K with which 
it is necessary to work in order to obtain the final solution. 

§ 13. The Equilibrium of Two Ellipsoids joined by a Weightless Pipe. 

In § 3 the problem is considered of the equilibrium of two masses of liquid, each 
constrainedly spherical, when joined by a pipe without weight. It was shown that 

VOL. CCVI. A. 2 E 
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the condition determining the ratio of the masses for given radius vector is expressed 
by a certain equation which was written /(a/r) = 0. Further, it was proved that if 
/(a/r) is positive the two spheres of Hquid are too far apart to admit of junction, and 
if it is negative they are too near. Finally we found that all these solutions were 
unstable. 

The solutions for the two spheres showed them to be always very close together, and 
as all the solutions for two ellipsoids, when they are in limiting stability, made them 
much further apart than were the two spheres, it seemed somewhat improbable that 
two ellipsoids could be similarly joined by a pipe, and cei-tain that they would be 
unstable if such junction were possible. Nevertheless, it seemed conceivable that the 
additional terms, which must appear in /(a/r) when the constraint to spherical form 
is removed, might alter the conditions so that the junction of ellipsoids by a pipe 
should become possible. It thus became expedient to solve a problem analogous to 
that of § 3 when the two masses of liquid are ellipsoidal. 

The conditions of equilibrium of two ellipsoids unjoined by a pipe are given in § 10, 
and the additional condition corresponding to junction by a pipe is 



dV 2. 2^1 
d\ ^ dX 



+ loy'~ = 0. 



In the present investigation I shall neglect the higher ellipticities, denoted // and 
jP/, and terms of higher order than those in 1/r^ 
With this degree of approximation we have 



V = (|upa7 



\' 



L(i+x)v"^"^"(i+x)^'^"^^^°(i+x)^^j 



+ {eE)2, 



where {eE)2 is given in (14) (with omission of terms in l/r'), 
Also 



"^ = 1 



00 



dv 



vo 



{v'--lf'{v'-l/Kj^'' 



and ^ has a symmetrical form in K^ and K. 

a^ 
We have besides co^ = f Trp -3 (1 + ^), where C is given in (24) ; and / is given in (22). 

The differentiation with respect to X and the subsequent re-arrangement of the 
equation are rather tedious, and I will not give the details of the operations. It may 
however be well to note that k, K^ are functions of X, and that 

# ^ 1 . c^^ ^ 1 ,p 

d\ 3X(1 + X)'^' d\ 3(1 + X) " 
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I find finally that the equation of condition is y ( - j = 0, where 






.2 



{i+\y"{i-\"') 



^{48-20X-(l9-15X)cos'y-(9-25X)cos'/8} 



+ |^{20-48X-(l5-19X)cos^r-(25-9X)cos'B} 

+ ^^{288-112X-(260-140X)cos'y-(204-196X)cos'/8 

+ (87-63X)cos*7+ (59-91X)cos*/8 

+ (30-70X)cos'/8cos^y} 

3aa* 
5 6 Or 

+ (63-87X) cos" r+ (91-59X) cos* B 

+ (70-30X)cos'Bcos'r} 



+ ||^g{112-288X-(140-260X)cos'r-(196-204X)cos'B 



+ ^^^ {40 (1-X) - (18-22X) cos^y- (22-18X) cos^ T 

- (14-26X) cos')8- (26-14\) cos' B 
+ 15 (1 ~\) (cos' y cos' r+ cos' /3 cos' B) 

+ (7-3\)cos'7Cos'B + (3-7X)cos'rcos'y8} (44). 

In this expression c and C are respectively the longest semi-axes of the two 
ellipsoids, which are pointed at one another. 

We may derive xjj and ^ from Legendre's tables of elliptic integrals for 

or we may expand the integrals in powers of /c" tan' y and obtain the approximate 
formula 

,^ = H[l + (fi-l)(l+iK'^+-,V*+/3V*...)-iAc'HanV(l+A'c'^+i^K'\..) 

+^2'^'* tan* y(l+-^K'^..)-9Vnan«y(l + ...)] . (45), 

where fl = —. — log'c • The formula for ^ is of course symmetrical. 

sm y cos y 

It should be noted that when the two ellipsoids reduce to spheres, we have 

2 E 2 
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Thus 



/(i. 



9. 

2 



1 + X^/^+X^ 



2/3 



(l+Xy/»(l + XV'')r 






This is the form obtained in the solution of the restricted problem of § 3. 

We conclude that if /(a/r), as expressed in (44), with values derived from any 
solution of the problem of the equilibrium of two ellipsoids unconnected by a pipe, 
is positive, the two figures are too far apart to admit of junction, and vice versd. I 
have in fact always found it positive, although always diminishing as r diminishes, so 
that junction would seem to be always impossible, at least so long as the approxima- 
tion retains any validity. This might indicate that there is no figure of equilibrium 
shaped like an hour-glass with a thin neck. However, I return to this subject in 
discussing numerical solutions, and in the summary of results. 

In the case where the two masses are equal, X=l, and the above formula for y*(a/r) 
fails by becoming indeterminate. As this is a case of especial interest, it must be 
considered. 

Since the two shapes are now exactly alike, we may take k, y, /3, k to define either 
of them. 

27 a 



When \ = 1, the first term of /(a/r) becomes 
The second term becomes 



4.2^/V 



3t//a(l + X^/''^+X^/^) ^ 9i//a 
2{l-h\y'{l-\'''') 2(l + X)^/^(l-f-X^/^) 4.2^/'^ 



3ai//(X-l) 



All the terms in 1/r^ are one of the two forms 



F 



c 



XC 



1-X'^' 



or 



G 






Now, 



& 



a 



X2/3 



(cos/3cos7f^(l+Xf^' 



C 



a 



(cosBcosr)'^'(l+X) 



2/3 



In the limit /3 = B, y = T, and X = I ; and we find that the first of these forms 
becomes ^Fc^, and the second —■ |(tc^. 

Again, of the terms in l/r\ those in c^ and (7^ are of one of the two forms 



F 



fc'-xcn 


y^ 


\\c'-0'l 


Li-x'^'J 


or G 





In the limit the first of these reduces to —^Fc*, and the second to — f (rc*. 
The last term in l/r^ has a common factor 1~X, when y = r, y8 = B, and 



1-X 



l+X'^'+X'" 



3 



{ 1 + Xf ( 1 - X^/8) ( 1 + X)i''^ ( 1 + X"'') 2. 2'" 
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By means of these transformations we find for \ = 1 



a 



f-) = 



,r 



16^/^ 



|-+'/'a+ ;^{-74 + 47cosV+67cos^;8) 



ac 



+ ^^(-272 + 300 cos^y + 356 cos' j8-123cos*y-151 cos^^S 



7 Or 



— 1 1 cos' yS cos' y) 



. (46). 



a 



If we put /3 = y = 0, and note that c becomes ~j^ , this expression reduces as 
before to the correct form. 



§ 14. Ellipsoidal Harmonic Deformations of the Third Order, 

The two ellipsoids whose forms have been determined are subject to further 

deformation by harmonic inequalities. The expression for the ellipticity f' 

corresponding to all the cosine-harmonic functions for i greater than 2 and s even 

is given in (25), viz. :— 

_ _ Jg/{r/A;}_t/{vo) i/ (1) (H:/ (^tt) ; 



//=:-H2^+l)- 



XE/^/-A/ 



k^ cos' y 



Z\.i^ sin'/3 _ 



1| 3(3 + .') F ■ 

14k' '• 



r 



■ (47). 



I shall begin by considering the ellipticities j's and fi corresponding to i 

5 = 0, 2. 

I define 

\%^{v) = v{K'v-q% (s = 0, 2) 



3 and 



(48): 



where q^ — f [l + /<:^+'v/(l — |/c^ + /c*)], with upper sign for 5 = and lower for 5 = 2. 

Then 

dv 



fKJJ 
f 
7/ 



[f / (v)]' (^'- ir (^'- 1 AT ' 

Since v is always greater than unity, the function under the integral sign may be 
expanded in powers of 1/v, as in (7), (8), (9) of § 5, and the integration may then be 
effected. In this way I find 



©/ (p) 



7k V 



2. .4 



1 + 



7 (1 + ff') + 10^' , 7[9(3/c'+2K'+3)+28g'(/c'+l)-f40(?^] 



2.9/c'r' 



+ 



8.9.11ffV* 



+ ... 



(49). 



In all the cases we have to consider k^ is nearly unity and k^ is small. Then, 
since g^^ = ^ [4 — 2/c^^+y/(l — /c^^ + 4fc^*)], and since the function under the square root 
may be expanded in powers of /c^^, we may obtain approximate expressions for (f in 
the two cases 5 = 0, 5 = 2. 

When 8 = 0, we have 
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When s = 2, 



9: 



1 



2/C -t-g'^ -T xqK , 



If we substitute these values for q^ in (49), and express the functions of k^ therein 
in terms of k^, we obtain the following results : — 



k 



an- _ 






fjr2~l] ^^ 3\^^ 2^ +32'^ +64'^ --V^ +33 V^ + '^ +28'^ ^^14^^ •••/ 4--- 



Shorter forms may be given to these by making the expansions run in powers of 
] Jkt^ ; we then have 



aa (| 



^\l+^H^A-Qx''-h-'*-\>c":.) 



¥ 



Kr 



+ 



¥ 



3 5 

3¥~0"* 
KT 



€: 






/^^ 



yi; 



2 



,2.^4 



/kV 






.•) 



^^ 



>- 



kV^ 



It will however suffice for our purposes to take 



V ■ 



(50). 



€3('^ 






y^^ 



7kV^ 



1 I 10 '^ 



1 + 



KT^ \ 



4 ^ 
3 ^^2 



1> 



(51). 



-/' 



The next task is to determine the product ^3* (^0) l^z (1) C3* (^tt) {s = 0, 2). 

The form of the ellipsoid is determined by z^o? where vq = — ; — = — — -z. . 

fcsmy smp 

If we write Ai^ = 1— g^^sin^y, with the definition of ^3* given above in (48), we 

have 

A/ 



^/ (^0) 



• • 



K sm" y 

It will be remembered that q has a different value according as 5 = or 2. 
I now make the following definitions, 



(52). 



so that 



P/ (/.) = II (Ky-q'), 0/ ((^) = {q''-K" cos^ <P) y(l -k'' cos^ <^) ; 

W(l)C/(j7r) = ^'^(K^-g^) (53). 



It is easy to show that rigorously 

q"{^-q') = -i-gll-K"-K"±{l-iK")^{l-K"+iK")} {s=0,2). 
Whence approximately, with the upper sign for s = 0, 

q'^^^-q^) = -ig(l-^'^+-igK'^+OK"...) (54). 
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With the lower sign for s = 2, 

g'^(K^-g^)=-lK'*(l+OK'^--T^/c'*...) (55). 

The three last expressions (53), (54), (55), give the two values of ^3'(1) Cs^d-Tr). 
We now turn to the functions ^Ds". The general definition of § 5 was 



®/ = ^^^ j [i/(/^) C/(<^)]> da, 



where M was the mass of the eUipsoid. Hence in the cases under consideration 



E; 



3 



7 sin^ ^ 



Arrk^ cos ^ cos y 



[W{l^)Os'{<]>)Jpd<r (s = 0,2). 



These integrals are evaluated in (38) of my paper on ' Integrals ' ; whence I find 

E/ = |i>[(l-i«'^)(l-'c'^-f'c'*) + (l-K'^ + f/c'*)i>], 



where D = ± y^(l— ac^^+4ac'^), with upper sign for s == and lower for s = 2. 

It will, however, suffice if we use the development of D in powers of ac'^ The 
result is, in fact, given in the equations next below (38) in * Integrals,' and they are 



Z, = (IY{1-2k''+^^K..), 



er/=3^'c'Ml-K'^ + fiK^..). 



(56). 



From (53), (54), (55), and (56) we now find 






_ 2 5 ^ K -f iq/C ... __ 



>v 



1-2k"+Uk'\.. 



W{^)Ci{hr)_ x^A + Ok" 



9^/4 



16 



®3' 



1-k"+Uk'\.. 



= ¥(i+'c'^-i'C^..), 



= -^(l + /c'^-f/c'^..). 



> 



• (57). 



J 



Thus from (51), (52), (57), we find 



r 



3\% 



ms(^m^^o)%{i)c^{i7T) = - 



25Ai' 



12Xsin''^ 



(l + /c'^-J/f'*...)^Yl + io^ 



2\ ^ 



9 ^1' 



> 



€/(T)i/(-o)i/(l)0/(i-) = ix||J^(l + -'^-t-'^)$(l + tS)- 



By (47) fi is equal to the above expressions divided by 



(58). 



J 



-A,'- 



F cos^ y 



1 + 



3(3 + K^)^ 



SXr'sin^^L 14k^ ^ 



...J. 



It therefore remains to determine ^3* and Ai\ 
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In order to find Ai^ we expand A in powers of k'^ tan^ y. Thus 

whence 

A/ = -^i-{[H-fK'^+i|K'1[l-cosV^]~i'^"tan^r(l+l'^") + i'f"fcan*r}- 
z sm y 

The result may, of course, be obtained with higher accuracy if it be desired to do so. 
By (25) of the ' Pear-shaped Figure ' 

sm y J Ai A 

where A/ = 1 — ^^ sin^ y , with appropriate values of g' for 5 = and 2 respectively. 

It is clear that under the integral sign 1/A may be expanded in powers of k'^ tan^y, 
but this is not possible with l/A/. 

The two cases 5 = 0, 2 have to be treated by different methods, and I begin with 
3, where 5 = 0. 

Since 

— [l + ^/c'Han'y + ffc^nan^y...], 



A cos y 
it is necessary to consider the integrals with respect to y of the functions 



sin^y sin^y sin^^y 



cos y A/ ' cos^ y Ai^ ' cos^ y Aj^ * 
Writing x = sin y, and q^ = — ^, this is seen to be equivalent to finding 



)(!-; 



a 



(1 -x') {a'-xy ' J (1 -xj {a'-x'f ' J (1 -x^)» {a'-xj ' 
Now 



{l-x'){cf-xy 2{a'-iy\l-x l+xj 4: {a^~l)\{a-xy {a+xf/ 



aU3a'-5)/ 1 , 1 



4c{a^—iy Xa—x a + 



X, 



^ i_ .^^:--3 /j_ , M+_i / i_^ , 1 



{l-xj{a'-xy ' 4 (a^-lf \l-x \+xl 4 (a^- 1)^(1 -x)^ {l+x)' 

a^3a''-7)/ 1 , 1 \ «« /I , 1 

~ ^-=^1 1 : -f , / ■>. 7T5 7 Vo + 



x"-' 



i{a?—lf\a—x a + xj i{a^—iy\{a—xf (a + x) 
- 1 I 3(21a^-18a'^+5) / 1 1 \ l7a'-9 ( 1 1 



{l-xy{aj'-xy 16(a^-l)* \l-x l+xj lQ{a''-lf\{l-xf {1+x) 



%{a^-lf\{l-xf [l + xf] 4.{a^-lf\a-x a+x, 

a^ / 1 1 \ 

4.{a'-lf\{a-xY {a+xfj ' 
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Each term in these expressions is integrable, and the limits of integration are 
X = siny to 0. After integration a mnst be put equal to 1/g, and since the result 
will only be true as far as k^^ we must put for q^ its approximate value f (1 — |-fc^^— f/c'^). 

When these processes are carried out, and the formulae combined, we find 



_ ^1 



4 



^ 4 ^* .5 



q^ sur y 



-i+f(i+iK'^+f|K'*)n-|(i+o^'^-3V*)-l 

--iV'secV(l+fM+-r2V'secV| • • « (59). 



In the coefficient and in Ai^ (or 1 —if sin^ y) it is as easy to use the rigorous value 
of (f as the approximate one, so it may be well to repeat that 

q^ = I [1 + k^-~v^(1-|ac^+/c^)], and n = -J- log,l±^Jt^. 

sm y cos y 

If a process parallel to that adopted for finding ^Is were adopted in the case of ^3^, 
it would lead to a divergent series, but fortunately a much simpler process is 
available. 

In the case of ^ = 3, 5 = 2, we have q^ = 1— |-/c'^-f f/c'^ ... , so that 

A/ = cos^ y [1 4-ifc'^ (1 -f k'^) tan^ y]. 

Accordingly we are now able to expand 1/Ai* in powers of k^^ tan^ y. 
We have then 

^ ^ -[l-K'^(l--fK'^)tan^y + |/c'HanV...], 



A cos y 
Whence 



^/ cos^ y 

^ ^ • [I-^k' tan^ y + 1 /c^^ tan^ y . . . ]. 



,:^6 ^. ^'6 



AAi^ cos^y^ 2 \ 2 / / 8 /J 



Now 






ol^^^== '^">'t« *^^'^-^*^"''>'+^*^^^*^-^'s^*^"'''^+tM^ ^''-')^' 



where 

n - _L loo- ^ + ^^^ y 

sm y cos y 
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Combining these formulse, I find 

^/ = ^4- [{tan* y-l tan^ y+\^ (O-l)} {1 + A^'^+-^3^V*} 

-'^''(l + 3V)tan^y+WtanV] . . (60). 

The form is somewhat awkward, but I have not been able to reduce it to any 
better shape. 

1 15. The Values of %i for Higher Harmonic Terms, 

For higher harmonic terms it is necessary to adopt the approximate forms of the 
functions investigated in ''Harmonics." The development is there carried out in 
powers of a parameter ^, which I will now write ^o to avoid confusion; this 

parameter is equal to ^ ? ^^ *^ ^//^+|■K'^.. of the present paper. 

The functions are here defined by 




but in the notation of §10 of "Harmonics" this would be called ®/(v)/E/. Thus, 
if [(!Ei' {vy\ denotes that function as defined in " Harmonics," we have 

@/ (v) = L'^^^V /J . 
We have for the approximate expression for p/ 

The investigation on p. 500 of "Harmonics" shows that the leading term of 



2H\ i+sl 



i + s-2] , o i+s + 2 ! 



~> 2i?T! v^ [1 +^0'?-^ "T^iTr "^^'•^'■^^ ^m 

This has to be divided by IE/, the formula for which is given in § 10 of " Harmonics," 
and we thus obtain the leading term of ©/ (v). 

For the second term it will suffice if we take /3o as zero, so that it is only necessary 

df 
to consider Qi {v), which is equal to {v^—iy^-r-^ Qi {v). 
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Since 



Qi (v) 



2' {i ly 

Z/h ~| J- 1 



1 ^ ^■+2! 

+ 



i+l 



2.1!^!(2i+3y 



+ 3 



"T . • . 



by differentiation, and by expansion of {y^—XY'^ in powers of \\v^ we obtain 






1 + 



(^'+2)(^^+l)+s' 
2 (2^ + 3) I/" 



+ ... 



Accordingly, .in order to find the second term to the degree of approximation 
adopted, it is merely necessary to multiply the leading term by 

1 + (^'+2)(t'+l)4-s ^ 
2(2^+3)1/^ "■ 

In order to find the leading term explicitly we have to insert for the g-'s their 
values, and after some tedious reductions I find 



<®, 



I, 



^^^^,|S{l+iA(^ + 2i) + o\A^[-2T+^^(.^+3) 



+ 2S(6^•-l + 2.-) + 2^•(3^•+l)-.^} x {l+ ^^+||^^+^3^ + ^^ gj • ■ (61), 



where 



S = ^' (^' + 1) T = (^-l)^(^+l)(^ + 2) 



This formula fails for the cases of s = and s = 2, and these cases have to be 
treated apart. Following a parallel procedure I find 



i'T 



^h 



^^^^^'t^ {] +iA(S + 20+d:^^„^[29S^+(l44*^ 



X \ 



1 + (i±ll(i±i)±f 1' 

2 (2*+ 3) V 



®^(0=|^J^{l-W(^-l) + li8A^^(^-l)(7*^-3^-6)}. 



2 1 



1 + (i±2)ii+i) ^ 

2(2^ + 3) r'J 



(62). 



(63). 



The values for i less than 3 are not required, and when ^ = 3 these formula are 
found to agree mutatis mutandis with those of the last section. 

It is pretty clear from general considerations that the higher inequalities 
corresponding to harmonics other than the zonal ones must be very small. I 
have, in fact, computed the third tesseral harmonic inequality {i = 3, 5 = 2), and 
find that it is so very minute compared with the third zonal inequality (i = 3, 5 = 0) 
as to be negligible. Accordingly it appeared to be a waste of time to develop 
formulae for any other than zonal inequalities for values of i greater than 3. Thus 
of the formulae just determined the only one of which actual use is made is (63). 

2 F 2 
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§ 16. The Fourth Zonal Harmonic Inequality. 
In developing the expressions for the higher harmonic inequalities it seems to 

1 2 

be most convenient to retain the parameter ^q, which is equal to ~ h, instead of 

1+/C 

developing in powers of k^^ as heretofore. 
On putting ^ = 4 in (63), we find 

^*\I/ ~ 57f:9;^{ l+TT^j (l"~3po + -i-6-po ). 

With the notation of ^^ Harmonics " we have 

C,(<^) = l-5;8o COS 2(^ + f|y8o' COS 4:<j>. 

Accordingly 

f ,(1) = P,(l) = 1, 

Again, from § 22 of " Harmonics" for type EEC, i = A, s = 0, we have 



6 



But Wi is this integral multiplied by 9 and divided by 3 times the volume of th 
ellipsoid. 
Therefore 

In this formula the coefficients of the powers of ^o increase with great rapidity, 
and the approximation may not be very satisfactory ; nevertheless it is the best 
attainable without an enormous increase of labour. 

Combining our several results 

It remains to find f 4(^0), and the denominator in the expression for /^ which 
involves ^4^—Ai^. 
We have 

dv 



= [f 4W] 



00 
2 



[f.(-)T(-^-ir(-^-i/'<r 



1 

Inside the integral sign I write v = — ; — - and change the independent variable 

^ ^ K sm \|/ 
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to xfj; I also put sin x = ^ ^^^ ^* ^^ the surface of the eUipsoid we have ^ = y? 

V = jS, and since pq = — = -- — 5? these are the values to be used in 1^4,{vo) outside 

K Sin y sm p 

the integral sign. 

Then we have 

dv 



f» 00 7 l*y 



= K 



sec X d^- 





Since 

P/(.) = i^(^^_i)[7(v^-i) + 6] = ^i{6 cot^x + 7 cot^xX 
P/(j.) = 105(t''-l)' =105cot*x- 



Since further 
we find 
where 



%{^) = P.(^)+iAP/(-)+ii8/3o^A*(^). 

P^(i/) = 1 + cs cot^ X + ^ cot* X, 



a=5(l+-|A), 

6 = ^/(l + 3A+i%/So^). 



^ . 



(65) 



y 



It must be noted further that when v = vq at the surface of the ellipsoid, ^ = ^8. 
It follows then that ^4 (^0) = l-\-co cot^ ^+& cot* jS ; and from (64) 



SX'^^^^\k r^'^'"' -^'^^^ -'^■^■■' 105Xrn ^^r 



( 1 - 8/80 + ^I^A') ( 1 + « cot' y8 + 6 cot* fi). 
It remains to consider the evaluation of ^4, which now assumes the form 

P / l+acot'/J + 6cot'^ Y sec Y rfii; 

It would no doubt be possible to split the subject of integration into partial 
fractions, and thus obtain an accurate value as was done in the case of ^3, but it does 
not seem worth while to undertake so heavy a task, because a sufficiently exact value 
may be obtained by quadratures. 

The method is employed in §18, p. 294, in my paper on '^ Stability,'^ and may be 
explained very shortly. 

I divide y into 10 or 12 equal parts — say 10 for brevity — -and let 8 = -j^y. 

I then compute eleven equidistant values of the subject of integration, say Uq^ Ui, 
... Uio, corresponding to ^ = 0, 28, 38, ... lOS. As a fact it is unnecessary to compute 
the first four of these, because they are practically zero. 

The equidistant values increase so rapidly that they are very inappropriate for the 
application of the rules of numerical quadratures. Accordingly I take an empirical 
and integrable function, say v, such that t'i2 = U12 and Vu — Un^ and apply the rules 
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of quadrature only to the differences u^—Vn. The result is a correction to the integral 
vdxjj, 

JO 

The empirical function which satisfies these conditions is 

tullog Hifl 
V = UiqC ' . 



log Hio 



When \]j = y = lOS, V = UiqI and when xjj :=: 98^ v — Ui^e ^''» = Uq, 
Then f . # = ,-,-1^ (1 -e-'.<-). 

In the cases we have to treat e"'^^^*'^^^''^"^"^^ is an extremely small fraction, so that 
practically \ vdxb = —--:: — ^7^- — j — r ; and this is the function to be corrected by the 

result of quadrature. 

For the quadratures we have 

^^10 = ^105 % = %, % = %o(~) . "i^i = y^wi — ) , &C. 

\%o/ K'^J 

Thus the equidistant values of the function to be integrated (arranged backwards) 
are 



0, 0, Us — Uio(-^] , Urj — Uioi-^] , &c. 



The first two are zero, the next three or four are sensible, and the rest are 
insensible ; thus the quadrature is very short. The correction is found to be very 
small, and we might perhaps have been content with the empirical integral without 
material loss of accuracy. 

§ 17. The Fifth Zonal Harmonic Inequality, 

This is treated exactly in the same way as the fourth, and I will only give the 
results. 



We have 



€5@.= Tro^TT^I(l+H^)(l-5/3o + ¥6W), 
C^{v) = v/(l-A cos 2<^).[l-7y8o COS 2<^+f|y8o' cos 4<^], 



Whence 



'' ^.fi)f.(l)0a(i7r) = - ^4^^!(l+ep)(l-13A + ^FA^)v/(l+^o). 



3X® 
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Then 



where 



Finally 



^^(x/o) = cosec/3(l + acot^/3+6cot^)8). 



^i 



"f 



cosec ^ (1 + g cot^ /3-^h cot^ ^)" 
cosec X (1 +^ cot^ ^+& cot"^ x)- 



sec X dy\f, 



which is to be evaluated by quadratures as was proposed for the fourth harmonic. 



§ 18. Moment of Momentum and Limiting Stability. 

The moment of momentum of the system is /w, but when we are determining the 
configuration of minimum moment of momentum, which is a figure of bifurcation and 
gives us the configuration of limiting stability, the conditions are different according 
as whether we are treating the problem of the figures of equilibrium where both 
masses are liquid, or Roche's problem in which the ellipsoid denoted by capital letters 
is rigid. 

Accordingly I write 



/ = f Trpa' 



1 

5 



\r 



TTx(^'+"") + (iTX?J^*"'''-' 



1 + X 



{B'^C% 



and for determining the angular momentum of the figures of equilibrium I take the 
whole expression for /, but for Koche's problem omit the last term. 



a 



Since M^ = 37^p-3(l + 0' I compute for Roche's problem 



r,3/2 



r 



+ 



\r 



L' 1 + X a' ■ (l+X)^a^J 



and for the figures of equilibrium 



a^/^ 



The moment of momentum is given by 



B'^C 



5 (1 + X)' a^ 



/o) = (|7r/)f'a^(/xi0r/x2). 

It will be observed that /xi and /xg are expressible by numbers for any given solution 

of the problem. 

Suppose now that we have a succession of solutions for equidistant values of y 
differing but little from one another. Then if the solutions lie close to the region of 
limiting stability, we shall find that one of them corresponds to minimum moment of 
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momentum, either of jjli or of 1^2, as the case may be. Such a solution is a figure of 
bifurcation, and of the two coalescent solutions one has one more degree of instability 
than the other. If one of the two is continuous with a stable solution, and if, 
moreover, in the passage to the undoubtedly stable solution it passes through no 
other point of bifurcation, one of our two solutions is secularly stable and the other 
unstable. 

Now, two liquid masses revolving about one another orbitally at an infinite distance 
are undoubtedly stable, and such a case is also continuous with one of our solutions. 
Further, Schwahzschild has proved that Roche's ellipsoid has no point of bifurcation 
from first to last, and as this is true of one such ellipsoid, it is true of two.^'" Hence 
we conclude that the minima of )Ui and [I2 will afford figures of limiting stability. 

§ 19. Approximate Solution of the Problem, 

It is clear that spherical harmonic analysis is applicable to the case when the two 
liquid masses are widely distant. When they are so much deformed by their 
interaction that that method becomes inapplicable, good results may be obtained from 
the formulae of the last sections by means of development in powers of sin y, and it is 
this plan which is especially considered in the present section. 

It appears from § 1 that when one of the masses is small compared with the other 
(X small), the configuration of limiting stability for the problem of figures of 
equilibrium occurs when the two masses are very far apart. As X increases, that 
configuration corresponds with diminishing radius vector. It seemed then probable 
that at least some of the solutions might be found by means of these series, and if 
this were so it might, in many cases, prove unnecessary to follow the same laborious 
procedure as in finding the limiting stability of Roche's ellipsoid. This view was 
found to be correct, and I therefore think it well to record the methods by which the 
developments may be obtained, without however giving the full details of the very 
laborious analysis. 

When the masses are far apart, the terms denoted e and rj in the equation for ac'^ 
and in that for a^r^ are small, and they must be neglected in the developments. 

Writing for brevity g = sin y, we may prove that 



0=.-^log,ii-^i^ = S-L^ tan^y- -l + tg'\ tan^ y = 1+2(^-1)^ 



sm y cos y 2n+l 



2n 







tan«y^-l + S(^-^)(^-^^-^% tan^ = ^^ ^ (n-l) (n--2) (n-3) ^,„^ ^^_ 

1.2 1.2.3 

Hence the developments may be obtained of the functions ctq, cfi, cr^ ... tq, ti, T2 ... , 

^ SCHWARZSCmLD, " Die Poiiicar^sche Theorie des Gleicligewiclits." 'Neue Annalen der k. Sternwarte 
Miinclien,' Band III, 1896. 
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and thence of vq, vi, vg... in series proceeding by powers of g^ ; and thence we may 
find k" in that form. 

The result as far as g* is 



k" = 



1 + X 
4 + X 



30 2 12(11 + 26X) 4 



7{l + X)\ 7'{i + xy 



I have also found the term in g*^, but shall make no use of it. 

With this value of /c'^ or of k^, which is 1 — k'^, we develop the expression for Si^/r 
in the same manner. The result is : — 



a' 



2 



3 



1 + X 

4 + X 



r 



1+ J + ^^^ 



88 + 53X 4 



7(4 + X)^' 7' (4 + X) 



By inversion we have 



2 . 2 5 (4 + X) a^ 

^ ^ 2(1+X)r^ 



25(5 + X)(4 + X) a^ 1 25 (69 + 2X)(2 + X) (4 + X) a^ 



2^.7 (1+xy 



v,6 



+ 



o- F COS iS cosy X 3 -i P n ii 4. F 

Snice r^~-7^^ — ^ = :: — r a' , it loUows that — 



2^.7^ (1 + Xy 

X \^/^ 



^.9 



a 



sin^/3 



1 + X 



sin^ /S \ 1 + X/ (cos ^ cos y )^' ' 



The semi-axes of the eUipsoid are given by 



¥ 



X \^/^ 



¥ c' 



a^ a^ sin^ y8 



1+X/ (cos /8 cos y) 



- = -2 cos' y, 3 = 3 



2/3 ' 2 



cos' ^. 



But cos'y = l—g^, cos'^ = 1— k'p', and therefore 



a 



2/3 



2 



1 + X. 



(l^g2>^-V3(l_.^2^2>^-l/3. 






1 + X. 



2/3 

) (1 



■cff'{i-KYY'^'; 



a^ 



/ X 



1 + X, 



f(i-/)-v^(i-Ky)n 



Setting apart the factor [X/(1+X)P^, which is common to all, these three are all 
expressible in the form, say 

i^= 1+ (^0 + ^1^')^^'+ {h-VhiK^^h^K^) g^ -^ (Co + Ci/c'+C2K^ + C3fC^)^^+ .e, , 

where the ao, a^^ h^, hi, &c., have different numerical values according to whichever of 
the three functions we are treating. 

Now the above formula for k^^ enables us to write 

K = ^0 + B^g + C^g . . , 5 

where the forms of .i4o, B^, Co are obvious. 
Hence we have 

jP = 1 + (^0 + ^1^0)5'^ + (&o + ^i-^o+?>i^o + Mo^)i7^ 

+ (Co + CiJo + C2^o^ + Mo^+&ijBo+2Mo5o + <^iC'o)5'^ .-. 
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In this way I find the following expressions for the semi-axes in series proceeding 
by powers of g^ or sin^ y : — 



c 



2/8 



1 + X, 



1 4" 



7 + X 



3(4 + X) 



g ™|- 



524 + 223X+14X^ 

3\7 {4+xy 



^ 



9' 






120926 + 86748X+ 19428X'+ 686X.^ 



3*.7'(4 + Xy 



9 



,6 






2/3 



1 ■\'Ki 



1- 



"T" -^A o 

9 



64 + oa+7a 4 

9 



3(4 + X)^ 3^7 (4 + X) 



?/ 



ajS 



a 



2 



\l+\ 



1- 



2 — A 2 

3(4 + X)^ 



_ 11122 + 2460X-1851XH196X =^ 

3*.7'(4 + Xf 

187 4" 11 OX"— 14a 4 
3^7 (4 +X)^ ^ 

28492 + 36723X + 14160X^-686X^ 



/. . . 



^ (66). 



3*.7'(4 + Xf 



9 •*• 



-/ 



It i^ easy to verify that the product of the three series is unity, as should be 
the case. 

The next step is to substitute for g^, g^^ o^ . . . their values in terms of a^^^^ and its 
powers. In this way I find 



c 
a^ 



2/3 



1 + X, 



1 JL, 5 (7 + X) a' 25(419 + 187X+11X") a 
6(H-X)r^ 2^.317 (1+Xf 



6 



^ 



^9 



+ 



125(998484-74769X+16503X^+488X')a' 



2^.3^7'{l+X)^ 



/p 



,9 



a" 



2/3 



a 



2 



1 + x. 



_ 5(5 + 2X) s? 25(11 + 37X-X^ ) a; 






9 



,9 



V 
a 



2/3 



2 



1 + X, 



6(1 + X) r^ 2^.317 (1 + X)^ 

_ 125 (23992 + 17895X+1587X'+178X') a 

2^3*.7'(1 + X) 

5 (2-X) a^ _ 25(157 + 119X-1 1X') a» 
6(1+X)r' 2^.3^7 (1 + X)^" r« 

125 (19114 + 23673X+11577X'-488X'') a' 

2^3*.7'(1 + Xf r 



} (s^'). 



y 



By writing 1/X for X we obtain the formulae for the axes of the other ellipsoid. 

The numerical coefficients increase rather rapidly so that the series are useless 
unless a/r is small, and accordingly this method fails to give any result for Roche's 
ellipsoid in limiting stability ; it is, however, useful for the problem of figures of 
equilibrium, as already stated. 

If we had relied on spherical harmonic analysis, we should only have obtained the 
terms in a^^^. 
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In order to obtain the expression for the angular momentum, which has to be a 
minimum for limiting stability, we must evaluate C Now, from (24) and (15), we 
have 

lOr 



( = -——,(2c^—a^—}f+ same in A, B, C) 

3 



+ 



56r^ 



[3 (a*+6*)+8c*-8c'(a'+6^) + 2a'62+ same in A, B, C\ 



+ ^[2(^V+5^6^+(7V) + (^^+5^+C^)(a^+?>^+c^)-5C'(a^+6^) 



2 Or' 



-5c^(^^+5^) + 5c^C^]. 



By means of the above series for a^, })\ (?, and of their analogues for A^, B^, C^, 
I find 



a 5 Q 8 «^ .10 



9 



.5 



/V»^ /V*-'-'-' 



where 



7 — ^» 



X^/^(7 + X) + 7X+l 



(1 + X) 



5/3 



m 



5 X^^''(82 + 38X+X^) + 82\"+38\+l 

"'*■ (1+XF 



w 



- 215 X^'^' (33 + 1QX+X") + 33X^+ 1QX+ 1 15 X'^^^(15 + 102X+15X") 

-224 (i+xyo/^ "•■ i« (l+X)i«/» 



(i+xy 



Now we have to evaluate the moment of momentum given above in § 18, viz., 



f.. = i|)"'(i+£r 



;ij^,(6'+o=)Hji-,(B"+c>)+ii:^^; 



When ¥, c^ J5^, (7^ have their values attributed to them, we find 



1^2 = 



a 



3/2 



where 



(1+X)^ \r/ 



(1+cy^^ 



1/2 



<ai" A^w rt'*'' 



^(1+X'^/^)(1 + X)v^ 



T = 



U = 



6(1 + X)^/''L 
5 

25 



X^/='(5 + 2X)+i(5X+2) 

X 



m • 

X'''^(131 + 34X+llX')+i(131X^+34X+ll) 

X 



X'/' (40367 + 25548X+ 2463X'+ 488X'') 



+ ^(40367X''+25548XH2463X+488) 
X 



In § 1, where the same problem is treated for two spheres, we had Z, m, n^ S, T, U, 
all zero. 

2 G 2 
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In order to find the minimum moment of momentmii for a given valne of X, I 
compute I, m, n, R, S, T, U, and assuming several equidistant values of r compute 

values of fi^ x ^—^-^ . When the coefficients are computed we very easily find the 

A 

value of r corresponding to the minimum. 

When that value of r is found, we are in a position to compute the axes of the 
two eUipsoids. 

For values of X less than ^ the results found in this way would be satisfactory, 
and for X = I they are, I think, adequate. Even for the case of X = 1 the result is 
not very remote from the truth, for whereas the correct result for the minimum of 
angular momentum is r/a = 2'638, the result derived from this approximate method 
is r/a = 2*51. But it would have been impossible to foresee that the result would be 
as good as it is. 

Part II.— Numerical Solutions. 
§ 20. Roche's Infinitesimal Ellipsoidal Satellite in Limiting Stability. 

We require to find the form of an infinitesimal satellite (so that X = 0) revolving 
in a circular orbit about a spherical planet. When this problem is solved we shall 
be able to see how far the solution will be affected when we allow the spherical 
planet to become oblate under the influence of a rotation of the same speed as that of 
the revolution of the infinitesimal satellite. This last is what I have called the 
modified form of Roche's problem. 

The planet being spherical and X being zei'o, the small terms 4 e, 'q vanish, so that 
our solution becomes rigorous. 

The angular momentum of the planet's axial rotation is to be omitted, and the 
satellite being infinitesimal the momentum of its axial rotation is zero. Thus the 
moment of momentum of the system varies as the square root of the satellite's radius 
vector, and minimum momentum coincides with minimum radius vector. 

The solution of the problem has been obtained in two ways : first by Legendre's 
tables of elliptic integrals, and secondly by means of the auxiliary tables given above. 
In the first method, I knew with fair approximation by various preliminary compu- 
tations the values of k and y which lay near to the required solution. Now there is 
a certain function of fc, y, say /(sin~^ k, y), which vanishes when the ellipsoid is a 
figure of equilibrium ; accordingly I computed by means of Legendre's tables the 
following eight values of /(sin""^ /c, y) for integral degrees of sin~^ k and y :■ — • 

f{77'', 57°) -= +0-0000878, /{rr"", 58°) = •-0-0000624 
/{78°, 59°) = +0-0000724, /{78°, 60°) = -0-0000785 
/(79% 61°) = +0-0000562, /(79°, 62°) = -0-0000939 
/(80°, 63°) = +0-0000408, /(80°, 64°) = -0-0001046 
(probably the last significant figure in each of these is inaccurate). 



FIGUEE AND STABILITY OF A LIQUID SATELLITE. 



229 



Interpolating from these we find four values satisfying / (sin ^ y, k) = 0, namely : — 
f(77°, 57°-5846) = 0, /(78°, 59°-4798) = 0, /(79°, 61=-3744) = 0, /(80°, 63°-2806)= 0. 
With these solutions I find 



sin ^ K 


'/■/a 


77° 


2-467860 


78° 


2-458191 


79° 


2-455446 


80° 


2-460289 



By formulae of interpolation the minimum of r occurs when sin"^ k = 78^*8756. Then 
by a second interpolation this value of k corresponds with y = 61°*13835 and the 
minimum value of r is 2*45539. We may take then y = 61"^ 8''3, k = sin 78"" 52^*5^ 
whence /3 = bd"" 14^*5. Since cosy = 0'4827, cos^ = 0*5114, the three axes of the 
ellipsoid are proportional to 10000, 5114, 4827; IIoohe gave the ratios 1000, 496, 
469, and the radius vector as 2'44, in place of 2'45539. 

Turning now to the second solution, I solved the problem by means of the auxiliary 
tables in two ways, namely, for y = 60°, 61°, 62° and also for y = 57°, 59°, 61°, 63°. 

They led to virtually identical results, viz., that the minimum of r is 2*45521, 
corresponding to y = 61° 8'*4, sin"^ k = 78° 52^*0. 

Finally the solution for Roche's limit and for the ratio of the axes of the ellipsoid 
in limiting stability may be taken to be as follows :-— 



7 

61° 8f 



sin ^ K 

78° 52' 



cos y 

0*4827 



cos^ 
0*5114 






with uncertainty of unity in the last place of decimals in r and of half a minute of 



arc in sin ^ k. 



We must next consider the modified form of Roche's problem, in which the large 
body or planet yields to centrifugal force and becomes an oblate ellipsoid of revolution. 
The approximate formulae of § 19 show that when X = co or when X = 0, 



1 = 



3. 

4 5 



m 



145 



^ — 2.2.A 
'«' — 2 2 4* 



Hence in this case 



3 
4 












The solution of the modified problem can only difier slightly from that just found 
when the planet is spherical, and therefore we may compute ^ with sufficient accuracy 
by means of the values of a/r already found. I accordingly computed I for 
y = 60°, 61°, 62°, and found that in each case I was very nearly equal to 0*0088. 

Now it is proved in the footnote to § 10 that, when X = and when the planet 
yields to centrifugal force, e = 77 = C; as the value of C is found with good approxi- 
mation, it is easy to compute r for these three values of y. I thus find that in the 
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modified problem, minimum radius vector, and therefore limiting stability, occurs 
when r = 2-457, y = 61° 12', «• = sin 78° 50', /S = 59° 17' ; the axes of the large body 



are determined by the approximate formulae of § 10 to be — = 



:2 = 1-0304, - = 0-9418. 

9j Qj 



a 



It appears then that the yielding of the planet to centrifugal force makes very little 
difference, as was to be expected. 

These results are included in the table given below of results for solutions of the 
modified problem of Roche with finite values of X. 

§ 21. Roche's Ellipsoidal Satellite^ of finite mass, in limiting stability, the planet 

being also ellipsoidal. 

This is the problem which I describe as the ^'modified" problem of Roche. It 
seemed unnecessary to carry out the computations for the smaller values of X, since 
they are sufiiciently represented by the case of the infinitesimal satellite where X is 
zero. I therefore begin with the case of X = 0*4 and pass on to X = 0*5, O'G, 0*7, 
0-8, 0-9, 1-0. 

It seems well to describe the process followed in one case as a type of all. It was, 
in general, possible either by extrapolation from neighbouring values of X, or by mere 
guessing, to begin with some values of {, e, rj and their co-relative functions E, H for 
the larger body, which were somewhere near the truth. With these we could compute 
r, fc, r, K with fair approximation ; thence values of ^, e, 77, E, H could be calculated 
with close accuracy and the computation could be repeated. It was, of course, a 
matter of conjecture as to what initial values of y would be found to embrace the 
region of minimum angular momentum. 

I will now describe the process for X = 0*4. Passing over the preliminary stages 
in which fairly good values were found, we begin with the following conjectural 
values : — 



sin~^ K 

r 

sin"-' K 
log r 



46^ 

68° 24'-8 

33° 13'-9 

50° 52'-0 

0-40245 



48°. 

69° 32'-0 

34° 18'-8 

51° 30'-l 

0-39594 



50°. 

70° 4l'-3 

85° 18'-5 

52° 7'-5 

0-39060 



r. 

sin"^ K 

r 

sin""^ K 
log r 



32°. 

50° 12'-0 

43° 50'-9 

67° 15'-0 

0-41071 



34°. 

51° 18'-8 

47° 24'-4 

69° ll'-8 

0-39775 



36°. 

52° 34'-9 

51° 34'-l 

71° 37'-0 

0-38726, 



whence I compute 



C 0-056259 0-064863 0-074219 
e 0-045435 0-048102 0-050256 
rj 0-32004 0-36426 0-40838 



^ 0-047883 0-062231 0-082010 
E 0-140179 0-174168 0-218540 
H 0-40539 0-52564 0-69272. 



By means of these and the auxiliary tables I find 

sin-i K 68° 24'-6 69° 31'-8 70° 40'-5 
logr 0-40240 0-39591 0-39046 



sm 



K 50° ll'-8 51° 18'-8 52° 33'-8 
logr 0-41068 0-39768 0-38693. 
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The computed values are so very close to the conjectural ones, in so far as they 
have been as yet computed, that we might be content, but in order to illustrate the 
process when the conjectures are less satisfactory, I proceed to the next stage. 

By far the greater part of the discrepancy between assumed and computed values 
(which in some cases was considerable) arises from error in the assumed values of r. 
Now assuming k and K to be correct, it is very easy to correct the results for a 
changed value of r. 

In this case I find 

corrected I 0*056274 0*064873 0*074271 corrected t, 0*047890 0*062253 0*082137 
e 0*04545 0*04811 0*05030 „ E 0*14020 0*17423 0*21889 

7? 0*32012 0-36431 0*40867 „ H 0*40546 0*52583 0*69385 



Recomputing 

sin""^ K unchanged sin"^ K unchanged 

corrected log r 0*40240 0*39591 0*39047 corrected log r 0*41069 0*39768 0*38695. 
By means of these we find two formulae of interpolation, namely : — 

^ = 2*4884-0*0342 l ^'^t^ ) + 0*0032 f^^^^ 
a \ 2 / \ 2 

^ = 2*4985-0*0685 (^^—^°\+0'0075 f^JlM'^^ 
a \ 2 / \ 2 

These two expressions may be equated to one another, and therefore we have the 
means of finding simultaneous values of y and F, and thence by another formula of 
interpolation of k and K. Hence I obtain 

y. 46°. 48°. 50°. F. 32°. 34°. 36°. 

F 33^ 13'*4 34^ l7'-9 35° ir-4 y 43° 47'-4 47° 25'*4 51° 34^*4 

sin"^ K 50° 52^*2 51° 29^*5 52° 6^*3 sin""^ k 67° 12'-1 69° 12'*3 71° 35^*6. 

Comparison with the initial values shows that the conjectures were very good. 

It now remains to compute the moment of momentum, and as we are dealing with 
Boche's problem the rotational momentum of the larger ellipsoid is not required. It 
follows that the values of F and K are not used, and since they are only required 
for finding the shape of the larger ellipsoid, there was no necessity for a high degree 
of accuracy in them. The moment of momentum is- represented by the quantity /xi 
of 5 18. I find then 

y. 46°. 48°. 50°. 

/xi 0*348640 0-348300 0*348519. 

By formulae of interpolation the minimum value occurs when y = 48°12'*9 and 

K = sin 69° 39'*1. The corresponding values are - = 2*4848, F = 34° 25'*6, 

a 
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K = sin 51° 33''5. The last step is to compute the axes of the two ellipsoids from 
the values of /c, y, K, P. 

Of course the numbers set out above make no claim to absolute accuracy, but the 
results tabulated below are, I believe, substantially correct. 

The unit of length employed is the radius of a sphere whose mass is equal to the 
mass of the whole system. If it were preferred to express the results in terms of 
the mean radius of the larger body, all linear results would have to be multiplied 
by {l + X)^/^. 

We may now collect the results in a tabular form, as follows : — 



SoLUTioisrs for Roche's Ellipsoid in Limiting Stability. 
The unit of length is the radius of a sphere whose mass is equal to the sum of the 



masses, ie., ahc+ABC = 1, and 



abc 
ABC 



= X. 



A. 


7- 


sin~i K. 


a. 


h. 


e. 


r. 


sm~i K. 


A. 


B. 


1 -030 


r. \ 





61° 12' 


78° 50' 


-482-^00 


0-511-4-00 


1-0 -^oo 






0-942 


1-030 


2-457 


0-4 


48° 13' 


69° 39' 


-562 


0-603 


0-843 


34° 25' 


51° 34' 


0-815 


0-886 


0-988 


2 '485* 


0-5 


46° 40' 


68° 12' 


0-597 


0-642 


0-870 


35° 59' 


54° 30' 


0- 


792 


0-860 


0-979 


2-484 


0-6 


45° 5' 


%&" 43' 


0-627 


0*674 


0-888 


37° 14' 


56° 41' 


0" 


772 


0-836 


0-969 


2-490 


0-7 


43° 38' 


65° 20' 


0-652 


0-701 


0-901 


3S° 9' 


68° 18' 


0- 


753 


-815 


0-958 


2-497 1 


0-8 


42° 26' 


64° 4' 


0-673 


0-725 


0-912 


38° 57' 


59° 39' 





■737 


0-796 


0-947 


2-502 


0-9 


41° 25' 


62° 58' 


0-691 


0-744 


0-921 


39° 40' 


60° 47' 





722 


0-778 


0-937 


2-508 


1-0 


40° 15' 


61° 43' 


0-708 


0-762 


0-927 


40° 15' 


61° 43' 


0-708 


0-762 


0-927 


2 -514 



The cases X = 0*4, 0*7, I'O are illustrated by figs. 2, 3, 4. The meaning of the 
dotted lines near the vertices of the smaller ellipsoid will be explained in the next 
section. 

The distance r—{c-{-C) is the interval between the vertices of the two ellipsoids; 
the following are the values, using, however, more places of decimals than are 
tabulated above :— 



A. 


r-{c,+ C). 


X. 


r-{e + C). 





1-030 


0-7 


0-638 


0-4 


0-653 


0-8 


0-643 


0-5 


0-635 


0-9 


0-650 


0-6 


0-633 


i-o 


0-660 



It is remarkable how very nearly constant the intervening space remains throughout 
a large range in the values of X. 

^ The values r = 2*485 for X = 0*4 and r = 2*484 for X = 0*5 represent 2*4848 and 2*4844 
respectively ; it is probable that the last significant figure in the former is a little too large and in the 
latter too small, and that it might have been more correct to invert the 2 * 485 and 2 * 484 in the table. 
I give the result, however, of the computation. 
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xiS= 886 




xb^mz 



X ^~ -Bee- 



Fig. 2. EoOHi's ellipsoid in limiting stability, when A = • 4. 



x/?--815 



x^=701 




X^=:701 



x^=-ai5 



Fig. 3. EOCHE^S ellipsoid in limiting stability, when X = * 7. 



X b^'im 



xB=7m 




x^ = 762 



X ^-=-762 



Fig. 4. Eociie's ellipsoids in limiting stability for equal masses. 
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§ 22. Harmonic Deformations of the Ellipsoids. 

When X is infinitely small, so that the liquid satellite is infinitely small, the 
harmonic deformations are evanescent, and the same is true when X is infinitely great. 
We saw in § 14 that there was reason to suppose that Roche's ellipsoid in limiting 
equilibrium might be more markedly deformed for values of X midway between zero 
and unity than for the latter value. I therefore determined the harmonic deforma- 
tions in the three cases X = 0*4, 07, I'O. 

The formulse for the ellipticities j^, fi^ f, f are given in §§ 14, 16, 17 and their 
values may be found for the ellipsoids in limiting stability tabulated in the last 
section. It appears that in every case the amount of deformation is small, and 
therefore it was sufficient to compute the normal deformation at the two extremities 
of the c semi-axis, that is to say, at the points nearest to and most remote from the 
other ellipsoid. At these points the normal displacement outwards may be denoted 
§c, with numbers affixed thereto so as to indicate to which harmonic it is due. 

The results may be given in a tabular form, but it may be well to remark that fi 
(second tesseral of the third harmonics) was only computed in two instances, because 
its effect was obviously quite negligible. For the same reason f^^ //, f^^ f^ were 
not computed at all. The following are the computed values of the ellipticities :• — 



\. 


04. 


0-7. 


1-0. 


A 


0-161 


0-234 


0-163 


fi 


0-127 


-0-102 


t 


/. 


0-0042 


0-0066 


0-0059 


f. 


0-0009 


0-0024 


0-0033 



We saw reason to suppose that the higher harmonics would be relatively more 
important for the larger values of X, and there is evidence of the general correctness 
of this view. 

The values of the ellipticities afford us no idea of the amount of the normal 
correction, and I therefore proceed to tabulate the values of 8c, the prolongations of 
the c semi-axis. 

Towards larger ellipsoid 

X, 0-4. 0-7. 1-0. 



Scg 


0-0193 


0-0283 


0-0194 


w 


0-00039 


0-00068 


? 

t 


8C4 


0-00472 


0-00886 


0-00911 


Scs 


0-00110 


0-00384 


0-00629 


:al 8c 


0-0255 


0-0417 


0-0348 


c 


0-8433 


0-9006 


0-9270 


c + 8c 


0-8688 


0-9423 


0-9618 


he 
c 


1/33 


1/22 


1/27 
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A. 


0-4. 


0-7. 


1-0. 


SCg 


-0-0 193 


0-0283 


-0-0194 


Sc/ 


-0-00039 


-0-00068 


2 

• 


Sci 


+ 0-00472 


+ 0-00886 


+ 0-00911 


Scg 


-0-00110 


-0-00384 


-0-00629 


Total Be 


-0-0161 


-0-0240 


-0-0166 


G 


0-8433 


0-9006 


0-9270 


G + Sc 


0-8272 


0-8766 


0-9104 


G 


1/52 


1/37 


1/56 



The last line in each division of this table has been given in order to show the 
relative importance of the total correction. It is clear that the ellipsoid remains a 
substantially correct solution. 

These corrections to the semi-major axes are indicated by dotted lines at the 
extremities of the longest axis of the smaller mass in figs. 2 and 3, and of both masses 
in fig. 4. 

We have in the last section tabulated r—(c+ (7), the distance between the two 
vertices. Now, although I have not calculated the deformations of the larger 
ellipsoid, it is pretty clear that they must bear to those of the smaller one approxi- 
mately the ratio of X to unity. Accepting this conjecture, we have for the 8(7 of the 
larger ellipsoid towards the smaller one the following values : — 



A. 



0-4. 



0'7. 



1-0. 



SO 0-010 0*029 0-035 

The distance between the two surfaces of liquid is clearly r—{c + Sc+C+SC), 

Thus we have 

A. 0'4. 0-7. 1-0. 



r— 



r~{c+C) 


0-653 


0-638 


0-660 


8g+W 


0-036 


0-071 


0-070 


(c+8c+(7+S(7) 


0-617 


0-567 


0-590 



§ 23. Certain Tests and Verifications, 

In order to test how nearly the solution for Si^/r^ by series in (40) of § 11 would 
agree with the solution (36) of § 10 in terms of the F elliptic integral, I computed for 
X = 0*7 the value of r/si in the two ways for three values of y, and found the 
following results :- 



X = 0-7. 






y. 42°. 


44°. 


46°. 


r/a by series 2-5355 


2-4888 


2-4500 


r/a by elliptic integrals 2-5359 


2-4881 


2-4495 


2 H 2 







236 SIR G. H. DARWIN ON THE 

The agreement seems to be as close as could be expected when five-figured logarithms 
are used. 

Certain terms in ^, as expressed in (24) of § 10, were neglected on the ground that 
they are of higher order than those retained. But it appears from the approximate 
solution in § 1 9 that the coefficients of the terms retained are themselves small, so 
that we are really only retaining terms of the same order as others which are 
neglected. 

The most important of the neglected terms in ^ is 

and this is a term of the seventh order. It seems, therefore, well to compute this in 
one case, and see how large a proportion it bears to the whole value. 
Since 

Using this value and the other approximate values given in §14, where /s is 
determined, I find that the neglected term is 

7Xsin^/3^^^^'' ^'^ VA ' kt' 



^ Ai F cos^y/ 3(3 + K^) F \ 
^' ^' SXr^sin^^V 14/c^ r^ ''V 

The numerical value of this, for the case of Roche's ellipsoid in limiting stability 
when X = 0*7, is found to be +0*0016. Now, the value of ^, as computed from the 
terms retained, was found to be 0*0677. Thus the neglected term is about one 42nd 
of the whole. The neglect then seems fairly justified. 

I thought it worth while to discover how far the modification of Roche's problem, 
whereby the larger body is ellipsoidal, afiects the result. I find that whereas it makes 
but little difference in the solution for any single assumed value of y, it does make 
a sensible difierence in the incidence of the minimum of angular momentum, and 
therefore of limiting stability. Thus, when X = 0*5, 1 found in one of my preliminary 
solutions for Roche's modified problem that limiting stability occurs when r/a — 2*49 
(the more correct value is 2*484), but when the larger body is a sphere it occurs 
when r/a = 2*35. Thus we see that ellipticity in the larger body induces instability 
at a greater distance than if it were spherical. This might have been conjectured 
from general considerations, 
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§ 24. Figures of Equilibrium of Two Masses in Limiting Stability, 

In this case both masses are liquid. We saw in § 2 that when one of the masses 
is infinitely small, stability only exists when the two are infinitely far apart. When 
X is less than 0*5 we may obtain fair results from the approximate investigation of 
§19, but for greater values of X it is necessary to employ the laborious method 
adopted in determining Roche's ellipsoid. 

When X ~ 0*5 I obtain the following approximate results for the two figures : — 



r = 


2-574 


a = 0-62, 


A = 0-81, 


b = 0-66, 


B = 0-87, 


G = 0-81, 


C= 0-95. 



It is probable that the value of r derived in this way is too small, 

For X = 0*4, I found r = 2'59, but did not calculate the axes. 

The only other case in which the problem has been solved is for equal masses, when 
X = 1. The two ellipsoids are exactly alike, and I find limiting stability occurs for 
the following values : — 



y=36M8', /c = sin59°33^ r = 2*638, a = 0*723, &= 0*771, c = 0*897; 



anci 



r-2o = 0*844, 



This is illustrated in fig. 5. 



xi?-771 



X 5=771 





x/?=771 xi9=-771 

Fig. 5. Two equal masses of liquid in limiting stability. 



§ 25. Unstable Figures of Equilibrium of Two Masses. 

When X = the figure of minimum radius vector, when the larger body is rigid, 
is also that of minimum angular momentum, but for larger values of X there is an 
ellipsoid considerably nearer to the larger body than that which possesses limiting 
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stability. I have only determined the ellipsoid of minimum radius vector in two 
cases, viz., when X = 0'8 and 1*0. 

When X = 0*8 I find minimum radius vector to be r = 2*36, whereas limiting 
stability occurs for r =2*50. When X^O'B, r=2*36, the ellipsoids are determined 
by the following data :— 

y = 5^ 20', K = sin Tl"" 51', T -= 46^ 10', K = sin 64° 20' ; 
whence 

a = 0*619, A = 0-705, 

h = 0-675, B = 0-774, 

c = 1-063, C= 1-018. 

When X = 1, the minimum radius vector occurs when y is about 54° and is then 
equal to 2-343, whereas limiting stability occurs when r=- 2*514. I have not 
computed the axes, since it suffices to learn that there is an ellipsoidal solution when 
the two masses are considerably nearer than is consistent with stabilitj^ 

As y increases, the ellipsoids get longer and longer, and it is interesting to inquire 
whether they increase in length with such rapidity that, notwithstanding the increase 
of r, the interval between the two vertices continues to decrease, or whether the 
increase of r annuls the simultaneous increase of c. 

The following table of values, computed with fair but not extreme accuracy, affords 

the answer to this question. 

X = 1. 

y, r. c. r-2c. Differences. 

44° 2-429 0-962 0*506 



46° 2-396 0-983 0*430 

48° 2-370 1-007 0-356 

50° 2-354 1-034 0-286 

52° 2-345 1-064 0*217 

54° 2-343 1-097 0149 

56° 2-350 1-184 0-082 

58'* 2-367 1176 0-015. 



76 

74 



70 



-69 



68 



67 



67 



The differences of r—2c hardly diminish at all, and it is clear that the next entry 
would be negative, or in other words the two figures would overlap. 

These results are obtained on the supposition that our approximation is adequate, 
but the small terms ^, €, 07, which are really infinite series, show signs of bad con- 
vergence as y increases, I think it probable that when we get to these extrenie 
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cases the convergence breaks down. It appears, however, justifiable to argue from 
these results that the unstable body continually elongates until its end coalesces with 
the other elongated body. I have no doubt but that the same holds true when the 
masses are unequal, and that we should always find r—(c~[-C) diminishing until the 
two meet. The poorness of the approximation of course would prevent us from 
making good drawings when coalescence is approaching. 

§ 26. On the Possibility of Joining the Two Masses by a Weightless Pipe. 

This subject is considered in § 13, and it is there shown that if a certain function 
written / (a/r), for a given solution of the figures of equilibrium of two detached 
masses of liquid, is positive, the two masses are too far apart to admit of equilibrium 
when joined by a pipe without weight — and conversely. 

Now I have computed / (a/r) in a number of cases of Roche's ellipsoids in limiting 
stable equilibrium, and have found it always to be decisively positive. 

The corresponding function for two spheres is given in (3) of § 3, and its first term 
is +a7'?"l When we compute it for two ellipsoids, we find the corresponding term to 
have become negative, and the additional terms, which are given in (44), § 13, are 
also negative. Hence / (a/r) is decidedly less for two ellipsoids than it is for two 
spheres of the same masses with the same radius vector. Thus the deformation of the 
two bodies tends in the direction of making it possible to join them by a pipe without 
weight, but it seems certain that in the cases of the Roche's ellipsoids in limiting 
stability such junction remains impossible. 

I also computed /(a/r) for the much elongated ellipsoids which are roughly 
computed in the last section and finally overlap, and always found /(a/r) to be 
positive, as far as the approximate formula went. The additional terms tend, 
however, more and more to cause / (a/r) to vanish, and the approximation becomes 
very imperfect. Now I believe, although I cannot prove it rigorously, that if we 
could obtain a more exact evaluation of the forms of these elongated ellipsoids, and if 
further a more exact value of /(a/r) were calculable, we should find /(a/r) vanishing 
near the stage when the computations would show the two ellipsoids to overlap. It 
therefore seems probable that there is a figure of equilibrium consisting of two 
elongated masses joined by a narrow neck. These ellipsoids are very unstable when 
detached, and, according to the principles of § 2, it seems inconceivable that junction 
by a neck of fluid could render them stable. 

Part IIL^ — Summary. 

Since the foregoing investigation may be read by mathematicians, while astronomers 
and physicists will perhaps wish to learn the nature of the conclusions arrived at, I 
shall devote this part of the paper to a general discussion of the subject, without 
reference to the mathematical processes used. 
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Two problems are solved here simultaneously ; for the analysis required for their 
solutions is almost identical^ although the principles involved are very distinct. 

We conceive that there are two detached masses of liquid in space which revolve 
about one another in a circular orbit without relative motion — -just as the moon 
revolves about the earth ; the determination of the shapes assumed by each mass, 
when in equilibrium, is common to both our problems. It is in the conditions which 
determine secular stability that the problem divides itself into two. 

One cause of instability in the system resides in the effect on each, body of the 
reaction on it of the frictionally resisted tides raised by it in the other. If now the 
larger of the two masses were rigid, while still possessing the same shape which it 
would have had if formed of liquid, the only effect on the orbital stability of the 
system would be due to the friction of the tides of the smaller mass generated by the 
attraction of the larger one. Investigation shows that in this case, as the two masses 
are bi'ought nearer and nearer together, instability would not supervene from tidal 
friction until the two masses were almost in contact ; but it is clear that the 
deformation of the figure of the liquid mass presents another possible cause of 
instability. In fact, instability, as due to the deformation of figure, will set in when 
the masses are still at a considerable distance apart. It amounts to exactly the same 
whether we consider the larger mass to be rigid, or whether we treat it as liquid and 
agree to disregard the instability which arises from the friction of the tides raised in 
it by the smaller body. Accordingly we may describe the stability just considered as 
'' partial," whilst full secular stability of both bodies will depend on the tidal friction 
of the larger mass also. 

The determination of the figure and partial stability of a liquid satellite (ie., apart 
from the effects of the tidal friction of the planet) is the problem of Koche. He, 
however, virtually regarded the planet as constrainedly a sphere, whilst in general I 
have treated it as an ellipsoid with the form of equilibrium. 

It has already been remarked that, as the radius vector of the satellite diminishes, 
partial instability first supervenes from the deformation of the smaller body. It 
therefore hardly seems worth while to consider the partial stability of a system in 
which the liquid satellite (hitherto described as the smaller body) is greater than the 
planet. We may merely remark that in this case the problem comes to differ very 
little from that involved in the determination of the full secular stability of two liquid 
masses ; for if we consider the case of a large liquid mass (the satellite) attended by a 
small body (the planet), it clearly makes very little difference in the result whether 
or not the tidal friction of the small body is included amongst the causes of instability. 

This being so, I have not thought it worth while to continue the solutions of 
EooHE^s problem (modified by allowing the planet to be deformed) to the cases in 
which the satellite is larger than the planet. The ratio of the masses of satellite to 
planet is denoted above by X, and the field examined by means of numerical solutions 
extends from X = to X = 1, while the part omitted extends from X = 1 to X = oo. 
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Tidal friction is a slowly acting cause of instability, and from the point of view of 
cosmical evolution the partial stability of Roche's e^ilipsoids is of eyen greater interest 
than the full secular stability of the system. 

The limiting stability of Roche's liquid satellite is determined by the consideration 
that the angular momentum of the system, exclusive of the rotational momentum of 
the planet, shall be a minimum. This exclusion of a portion of the momentum of the 
whole system corresponds with the fact that we are to disregard the tidal friction of 
the planet as a cause of instability. If all possible cases of the liquid satellite be 
arranged in order of the corresponding (partial) angular momentum of the system, it 
is clear that for given momentum there Avill in general be two forms of satellite ; but 
when the momentum is a minimum the two series coalesce. If then we proceed in 
order of increasing momentum, the configuration of minimum is the starting point of 
two series of figures ; it is a figure of bifurcation, and one of the two series has one 
fewer degrees of instability than the other. 

One of the two series is continuous with the case of a liquid satellite revolving 
orbitally at an infinite distance from its planet, and this is a stable configuration. 
Moreover, M. Schwarzschild has shown^' that the whole series of Roche's ellipsoids 
does not pass through any other form of bifurcation. Hence we conclude that of the 
two series which start from the configuration of minimum momentum, one is stable 
and the other unstable. 

The unstable series of solutions is continuous with a quasi-ellipsoidal satellite, 
infinitely elongated along the radius vector of the orbit, and the radius vector itself 
is infinite. Since two portions of matter cannot occupy the same space, the infinite 
elongation of the satellite would be physically impossible, unless the order of infinity 
of the radius vector were greater than that of the longest axis of the satellite. Now 
it appears from the numerical results of § 25 that this is not the case, and that the 
satellite becomes more rapidly elongated than the radius vector increases. Hence if 
the solution of the problem were exact we should reach a stage at which the two 
masses of liquid would overlap. I shall endeavour hereafter to consider the 
interpretation which should be put on this result. 

A series of solutions for Roche's ellipsoid in limiting stability is tabulated in § 21, 
and the table gives the radius vector and the three semi-axes of each body. The 
unit of length adopted is the radius of a sphere whose volume is equal to the sum of 
the volumes of the two masses. Three of these solutions are illustrated in figs. 2, 3, 4. 
The section shown is that passing through the axis of rotation and the two centres, 
but the places are marked which the extremities of the mean axes would reach if the 
section had been taken at right angles to the axis of rotation. 

The table of § 21 shows that the radius vector at which instability sets in only 
changes from 2*457 to 2*514, whilst X, the ratio of the mass of the satellite to that of 
the planet, changes from zero to unity. The distance between the vertices of the 

"^ See reference in § 1 8 above. 
T0L. CCVI. — A. 2 I 
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two ellipsoids also remains wonderfully nearly constant throughout a wide range of 
change in the value of X; for when \= 0*4 it is 0'653, and when X = 1 it is 0'660, 
only falling to 0*633 at its minimum. 

Thus far I have been speaking of the modified problem of Roche in which the 
planet assumes the appropriate figure of equilibrium, but I have also obtained the 
solution of Roche's problem for an infinitely small satellite and a spherical planet. 
As stated in the Preface, the radius vector of limiting stability, which has been called 
''Roche's limit," is found to be 2*4553, and the axes of the critical ellipsoid are 
proportional to the numbers 10000, 5114, 4827. These maybe compared with the 
2*44 and 1000, 496, 469 determined by Roche himself When we consider the 
methods which he employed, we must be struck with the closeness to accuracy to 
which, he attained. 

For the infinitely small satellite ''the modification" of Roche's problem hardly 
introduces any sensible change in the results, but for satellites of finite mass stability 
will continue to subsist for a slightly smaller radius vector for the spherical than for the 
ellipsoidal planet. In other words, the ellipticity of the planet induces instability 
earlier than would be otherwise the case. 

Roche did not attempt to investigate how closely his equations were capable of 
giving the ellipsoid most nearly representative of the truth, nor did he estimate how 
far the ellipsoid is an accurate solution. These points are considered above, and it 
was the desirability of making the investigation with a closer degree of accuracy which 
occasioned many of the difficulties encountered. 

For the infinitely small satellite the ellipsoidal solution is exact, and with a 
spherical planet, but not for an ellipsoidal one, Roche's equations give that ellipsoid 
exactly. In this case, however, the change introduced by the modification of Roche's 
problem is quite unimportant. 

For finite satellites Roche's equations require sensible modification, and the solution 
of the " modified" problem is different from that of the unmodified one, although not to 
an important extent. But the ellipsoid derived from the corrected equations is deformed 
by an infinite series of ellij)Soidal harmonic deformations, beginning with terms of the 
third order. Of these, the only ones which have any sensible effect are those which 
may be described as zonal with respect to the satellite's radius vector. 

By far the most important of these is the third zonal harmonic, whereby the 
satellite assumes a somewhat pear-shaped figure, being sharpened towards the stalk 
end of the pear pointing towards the planet, and bluntened at the other end. In 
consequence of this deformation the shape is slightly flattened between the stalk 
and the middle. 

The fifth and successive odd zonal harmonics accentuate the sharpening of the 
stalk and the bluntening of the remote end. The fourth, sixth, and successive even 
harmonics also accentuate the protrusion of the stalk, but tend to fill up the deficiency 
at the remote end. 
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The general effect must be very like what results from the second approxnuation 
to the pear-shaped figure of equilibrium,^ for I found that the ellipsoidal form was 
but slightly changed over the greater part of the periphery, whilst a protrusion 
occurred at one end — in this present case pointing towards the planet. 

In figs. 2, 3, 4, the protrusions at one end and the bluntening at the other, as 
computed from the third, fourth, and fifth harmonics, are indicated by dotted lines. 
It appears from these figures that, at least up to the point when instability sets in, 
the ellipsoid remains surprisingly near to the correct solution. 

For an infinitely small satellite minimum radius vector also gives minimum angul%r 
momentum, so that the closest possible satellite is also in a state of limiting stability. 
But this is not the case for finite satellites, and there exists an unstable ellipsoidal 
satellite with smaller radius vector than is consistent with stabihty. Thus for a 
satellite of four-fifths of the mass of the planet the minimum radius vector is 
2*36, whilst stability ceases at a distance of 2*50. Again, for equal masses 
stability ceases at 2*514, whilst the possibility of an ellipsoidal solution extends 
to 2-343. 

If we follow the forms of the more and more elongated satellites, when the radius 
vector has begun to increase again, we find explicitly in the case of equal masses, 
and with practicaJ certainty for all ratios of masses, that the distance between the 
two vertices continues to diminish and finally becomes negative. At this stage the 
two masses overlap, a conclusion which is, of course, physically impossible. But 
the calculation is based on the assumed adequacy of the approximations, and it is 
certain that the harmonic deformations of the ellipsoids increase rapidly, so that 
each body puts out a protrusion towards the other. The two masses of liquid must 
therefore really meet before we reach the stage of overlapping ellipsoids. As far as 
can be seen, the approximation has become very imperfect— perhaps evanescent — 
before the two ellipsoids cross. It will be best to continue the discussion of the 
meaning of this result after we have considered the true secular stability of the two 
masses of liquid. 

If a satellite, being a particle, revolves about a rotating planet, whose tides are 
subject to friction, there are, for given angular momentum, two configurations (if any) 
in which the planet always presents the same face to the satellite. In one of these, 
which is unstable, the satellite is close to the planet ; in the other, which is stable, it 
is remote, t If the angular momentum of the system be diminished, the radius vector 
of the sta^ble configuration diminishes and that of the unstable one increases until the 
two coalesce. For yet smaller angular momentum there is no configuration possible 
in which the planet shall always present the same face to the satellite. We see then 
that amongst all possible configurations in which the planet presents the same face to 

■^ See ' Stability,' referred to in the Preface. 

t See 'Roy. Soc. Proc./ No. 197, 1879, or Appendix Gr (b) to vol. II. of Thomson and Tafi's 'Natural 
Philosophy.' 

2 I 2 
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the satellite, that one is in limiting stability, in which the two solutions coalesce with 
minimum angular momentum. 

A rotating liquid planet will continue to repel its satellite so long as it has any 
rotational momentum to transfer to the orbital momentum of the satellite. Hence an 
infinitesimal satellite will be repelled to infinity, and the configuration of limiting 
stability for an infinitesimal satellite attending a planet, which always presents the 
same face to it, is one with infinite radius vector. 

Very nearly the same conditions hold good when both planet and satellite are 
subject to frictional tides. In § 2 it is proved that when each body is constrainedly 
spherical, the radius vector of limiting stability is infinite when the ratio of the 
masses is infinitely small. The radius vector decreases with great rapidity as the 
ratio of the masses increases, and when the masses are equal, the radius vector of 
limiting stability is 1*738 times the radius of a sphere whose mass is equal to the sum 
of the masses, or is 2*19 times the radius of either of the two spheres. Thus, when 
the ratio of the masses falls from zero to unity (and this embraces all possible cases), 
limiting stability occurs with a radius vector which falls from infinity until the two 
spheres are only just clear of one another. 

When we pass from the case of the two spheres to that of two masses, each of 
which is a figure of equilibrium under the attraction of itself and its companion, and 
subject to centrifugal force, the calculation becomes exceedingly complicated. Since 
the radius vector of limiting stability in every case must be greater than that of 
Roche's ellipsoid in limiting stability, and since in the latter case instability sets in 
through the deformation of the smaller body, it follows that in every case of true 
limiting secular stability of the system, instability supervenes through tidal friction. 

When the ratio of the masses is small, we have seen that limiting stability occurs 
when the two masses are far apart. In this case the deformations of figur^e are small, 
and could easily be computed by spherical harmonic analysis. 

For finite values of the ratio of masses, when spherical harmonic analysis would 
fail, a fair degree of exactness in the result may be obtained from the approximate 
formula of § 19. There would be no serious error from this formula when the ratio 
of masses is less than a half, but for greater values of the ratio it seems necessary to 
have recourse to the laborious processes employed in determining Egohe's ellipsoids. 
I thought, then, that it might suffice to compute the configuration of true secular 
limiting stability in the case of equal masses. It is illustrated in fig. 5, and we see 
that the radius vector is 2 '63 8. We found that for a pair of equal spheres, instability 
only set in when the radius vector, measured in the same unit, was 1738. Thus the 
deformations of the two masses forbid them to approach as near to one another as if 
they were spheres. It should be noted that instability in this case must arise from 
tidal friction, because Roche's ellipsoid in limiting stability was found to have a 
radius vector of 2*514. 

When Poincar:^ announced that there is a figure of equilibrium bearing some 



FIGURE AND STABILITY OF A LIQUID SATELLITE. 245 

resemblance to a pear, he also conjectured that the constriction between the stalk and 
the middle of the pear might become developed until it became a thin neck of liquid 
joining two bulbs, and that yet further the neck might break and the two masses 
become detached. References to my own papers on this pear-shaped figure and its 
stability are given above in the preface, and the present investigation was undertaken 
in the hope that a revision of Roche's work would throw some light on the figure 
when the constriction has developed into a thin neck of liquid. 

As a preliminary to greater exactness, I have in § 3 considered the motion of two 
masses of liquid, each constrainedly spherical, and joined to one another by a weight- 
less pipe. Through such a pipe liquid can pass from one sphere to the other, and it 
will continue to do so until, for given radius vector, the spheres bear some definite 
ratio to one another ; or, to state the matter otherwise, two spherical masses of given 
ratio, revolving in a circular orbit without relative motion, can be started with some 
definite radius vector so that liquid will not flow from one to the other. 

In this system the ratio of the masses and the radius vector are the only 
parameters, and I find that the condition of equilibrium is a cubic equation in 
the radius vector with coefficients which are functions of the ratio of the masses. 
The cubic has three i^eal roots of which only one has a physical meaning, and the 
solution is illustrated graphically in fig. 1. The single circle on the right is the 
larger sphere, and it is maintained of constant size for convenience of illustration. 
The smaller circles on the left represent the solutions for various ratios of masses, 
which are the cubes of the numbers written on the successive circles. 

The solution of this problem seems to me very curious, but it does not possess 
much physical interest, since it is proved in § 3 that all the solutions are unstable. 

The distance between the two masses is much smaller than is the case with any 
of Roche's ellipsoids, even with minimum radius vector, and accordingly it did not 
seem probable that the parallel problem, when the two masses are liquid and 
deformed, would possess any solution at all ; nevertheless, it was worth while to 
pursue the investigation to the end. 

When the masses are ellipsoidal and are joined by a weightless pipe, the solution 
would become very complicated, but the question may be attacked indirectly. 
When the masses are spherical there is a certain function of the radius vector and 
of the ratio of the masses which must vanish when a channel, of communication is 
opened between them. If this function be computed for two given spherical masses 
with given radius vector, we find that it is negative if the two masses are too 
close together to admit of junction by a pipe without disturbance of their relative 
masses, and that it is positive if they are too far apart. 

When the figures of equilibrium of two detached masses of liquid are determined, 
it is possible to form the corresponding function, but part of it consists of an infinite 
series of which it is only practically possible to give the first few terms. Now I 
have computed this function in a number of cases of Roche's ellipsoids, and have 
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found that the few terms of the infinite serie>s are small, that the series is apparently 
rapidly convergent, and that the function is decisively positive. We may conchide 
then that in none of the cases, for which numerical results have been given, is it even 
approximately possible to make a junction between the masses ; and even if we 
could do so, the system would be unstable, because removal of a constraint may 
destroy but cannot impart stability. To find any possible solution we must consider 
cases where the two masses are much closer together. 

I think, however, that there must be a figure of the kind sought, for the following 
reasons : If the function referred to above be formed for given radius vector and 
ratio of masses, we find that its value is very much less than if the two masses are 
spherical. Thus the tendency of liquid to flow from the larger to the smaller mass 
(when they are too far apart) is much less than if the two masses were spherical. 
Every increase of ellipticity in the ellipsoids tends to diminish the function, and the 
series tends to become less convergent ; and besides I have made no attempt to 
evaluate the terms in the function which correspond to the harmonic inequalities 
of the ellipsoids, and these would tend to diminish the function still further. 

It was remarked above that two much elongated ellipsoids seem to coalesce finally, 
but that the approximations were not satisfactory. I find, however, that even to the 
end the function, as far as it could be computed, was still positive although much 
diminished. It appears to me then probable that if we could obtain a miore complete 
expression for the function, we should find that it vanishes before the two ellipsoids 
overlap. There is then some reason to believe in the existence of a figure of 
equilibrium consisting of two quasi-ellipsoids joined by a narrow neck ; but such a 
figure must be unstable. 

I have, in fig. 6, made a highly conjectural drawing of such a figure where the two 




Fig. 6. Conjectural drawing of unstable figure of two equal masses of liquid just in contact. 



bulbs are equal. The data are derived from the computations for the much elongated 
ellipsoids just before they are found to overlap. 

Mr. Jeans has considered the equilibrium and stability of infinite rotating 
cylinders of liquid. This is the two-dimensional analogue of the three-dimensional 
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problem."^ He finds solutions perfectly analogous to MACLA.umN's and Jacobi's 
ellipsoids and to the pear-sbaped figure. In consequence of tbe greater simplicity of 
tbe conditions, he is able to follow the development of the cylinder of pear-shaped 
section until the neck joining the two parts has become quite thin. His analysis, 
besides, points to the rupture of the neck, although the method fails to afford the 
actual shapes and dimensions in this last stage of development. 

He is able to prove conclusively that the cylinder of pear-shaped section is stable, 
and it is important in connection with our present investigation to note that he finds 
no evidence of any break in the stability of that cylinder up to its division into two 
parts. 

The stability of Maclaurin's and of the shorter Jacobian ellipsoids is, of course, 
well established, and I imagined that the pear-shaped figure with incipient furrowing 
was also proved to be stable. But M. Liapounoff now states! that he is able to 
prove the pear-shaped figure to be unstable from the beginning, and he attributes the 
discrepancy between our conclusions to the fact that my result depended on the 
supposed rapid convergency of an infinite series, of which only a few terms were 
computed. The terms computed diminish rapidly, and it seemed to me evident that 
the rapid diminution must continue, so that I feel unable to accept the hypothesis 
that the sum of the neglected terms could possibly amount to the very considerable 
total which would be necessary to reverse my conclusion. I am, therefore, still of 
opinion that the pear-shaped figure is stable at the beginning ; and this view receives 
a powerful confirmation from Mr. Jeans's researches. The final decision must await 
the publication of M. Liapounoff's investigation. 

But there is another difficulty raised by the present paper. 1 had fully expected to 
find an approximation to a stable figure consisting of two bulbs joined by a thin neck, 
but while my work indicates the existence of such a' figure, it seems to me, at present, 
conclusive against its stability. The weightless pipe joining two bulbs of fluid is 
clearly only a crude representative of a neck of fluid, but I find it hard to imagine 
that it is so very imperfect that the reality should be stable, while the representation 
is unstable. My present investigation shows that two quasi-ellipsoids just detached 
from one another do not possess secular stability. The vertices of such bodies 
would be blunt points nearly in contact; the introduction of a short pipe without 
weight between these blunt points would differ exceedingly little from two sharp 
points actually in contact. Is it possible that the difference would produce all the 
change from great instability even to limiting stability ? The opening of a channel 
between the two masses is the removal of a constraint ; the system does not possess 
true secular stability when the channel is closed, and we should have to believe that 
the removal of a constraint induces stability ; and this is, I think, impossible. 

If, then, Mr. Jeans is right in believing in the stable transition from the single 

■^ "On the Equilibrium of Eotating Liquid Cylinders," 'Phil. Trans./ A, vol. 200, pp. 67-104. 
t See reference in Preface. 
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cylinder to two revolving about one another, and if I am correct now, the two 
problems must part company at some undetermined stage. M, Liapounoff will no 
doubt contend that it is at the beginning of the pear-shaped series, but for the 
present I should disagree with such an opinion. 

I have no suggestion to make as to the stage at which the pear-shaped figure may 
become unstable, or as to the figure which must be coalescent with it when instability 
supervenes. These points must await the elucidation which they will no doubt 
receive from future investigations. 

One question remains : If my present conclusions are correct, do they entirely 
destroy the apj^licability of this group of. ideas to the explanation of the birth of 
satellites or of double stars 1 I think not, for we see how a tendency to fission arises, 
and it is not impossible that a period of turbulence may naturally supervene in the 
process of separation. Finally, as Mr. Jeans points out, heterogeneity of density 
introduces new and important differences in the conditions. 
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